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ABSTRACT 



The conventional low-power CW doppler projectile tracking radar is investigated 
through simulation. The power requirement for full range tracking with real-time track 
data processing stopped within 160 m of the point of impact is obtained. This power 
can be lowered by a factor of more than 4 if a set of three switchable filters is utilized. 
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I. IjNTRODUCTION 



At a proving ground, it is always desirable to record the trajectories of test-fired 
munitions for later analysis. This is especially true when anomalies occur. Continuous 
wave (CW) radars have often been adopted for this purpose because of their simplicity. 
Compared to a pulsed radar of similar range performance, a CW radar has no modulator 
in its transmitter, hence no range gating is necessary in its receiver. Furthermore, a CW 
radar has a bandwidth corresponding to the doppler shift, which rarely exceeds 100 
kilo-Hertz (kHz) while a pulsed radar has a band width in the mega-Hertz range. Thus 
a CW radar can be designed to cany out coherent signal integration over a longer time 
interval. This increases the integration gain and lowers the required transmitter power 
for the CW radar. 

High gain antennas of ver}' narrow beamwidths are also used to reduce the trans- 
mitter power requirement. It is important to aim such antennas accurately on the target. 
To acquire the hearing and elevation angles of a target directly, several receiving feeds 
together with complicated comparator and multiple data channels are needed. Since the 
trajectory of a test round can be roughly predicted, an alternative to the actual acquisi- 
tion of target tracking angles is by steering the antenna according to the predicted tra- 
jectorx'. The directly measurable data will be the doppler shift of the target as a function 
of flight time. Discrepancies between this data and the doppler shift based on the pre- 
dicted trajectorx' should be carefully analyzed and explained. When an anomaly occurs 
which causes the target to break off track, this recorded data will be valuable for diag- 
nosing the problem. 

■Many CW radar tracking systems are available. For example, a model DR 5S2 X- 
band Instrumentation Radar manufactured by TERM.A Elektronik AS of Denmark is 
being installed at the Ta-Fu Proving Ground in Taiwan. These systems typically track 
a long-range projectile only up to the summit of the trajectory. Since the point of impact 
of a long-range projectile is of great interest, it is desirable to extend the range of the 
radar to cover as large a portion of the trajectory as possible while keeping the system 
relatively simple. The purposes of this thesis are to investigate the maximum range at- 
tainable using a simple CW radar and to determine the required minimum transmitter 
power. 
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When tracking a projectile, a CW radar detects only the doppler frequency shift due 
to the projectile motion. This frequency shift is proportional to the component of the 
projectile velocity pointing radially away from the radar. By combining this information 
with the tracking angles, the ground range of the projectile can be deduced through in- 
tegrating this velocity component. A low-pass filter can be employed to perform this 
signal integration in real time. The output can be used to update the predicted track if 
a data filtering mechanism such as the Kalman filter is utilized in real lime. Since in- 
formation on the doppler integrating filter is not provided by any manufacturer, one way 
to implement this low pass filter is suggested in Chapter 2. The design described therein 
is consistent with the operational instruction for the aforementioned system. It is needed 
for tracking simulation required in Chapter IV. The trajectory-updating Kalman filter 
is a desirable improvement to currently available systems on the market. Its design is 
outside the scope of this thesis. 

As e.xplained above, trajectory estimation is necessary for tracking with a simple CW 
radar. Since the radar and the gun are moving with the spinning Earth, the familiar 
reference frame in which the radar and the gun are stationary is an accelerating frame 
which is difficult to work with. Instead, by ignoring the rotation of the Earth around 
the sun, the inertial frame of reference having its origin at the center of the Earth and 
its z-axis along the spin axis of the Earth simplifies the equations of motion of the 
projectile. The z-x plane can be chosen to contain the radar at i = 0 when the projectile 
is fired. To obtain the familiar ground-based quantities, transformation to the inertial 
frame in which the radar (or the gun) is passing instantaneously through the z-x plane 
can be performed. In this frame, the distance and the relative velocity of the projectile 
from the radar (or the gun) can be found easily. 

The projectile will be spinning around its a.xis if a rifled gun barrel is used. The lift 
force due to air pressure will rotate the projectile to align its axis along its air velocity 
vector. The kinematics of the projectile, including its angular momentum, is most easily 
described in terms of the motion of, and the rotation around, its center of mass [Ref 
1]. The orientation of the projectile, which lies along one of its principal axis of inertia, 
can be described in terms of two of the three Euler angles [Ref 2], leaving the third angle 
for measuring the spin of the projectile. 

It can be shown that the effect of gravitational force on the projectile orientation is 
negligible. The ratio between the gravitational force on the center of mass to the force 
affecting the angular momentum is of the order of the square of the ratio of projectile 
size to the Earth radius [Ref 3|. Thus the change of the orientation of the projectile axis 



is due largely to the resistance of air. For a streamlined body, the lift from the air flow 
is large when there is an angle of attack [Ref. 4]. Furthermore, the rotation of the 
projectile axis is slow. It is reasonable to assume that the orientation of the axis of the 
projectile is always along its air velocity. This assumption will be utilized to determine 
the radar cross-section of the projectile. 

It is easily seen in the ground coordinate that the projectile velocity is much greater 
than the wind velocity so that the latter only introduces small corrections to the 
projectile axis. Since the wind acts on the projectile over the entire trajectorx'. a better 
estimate of the range can be obtained by including the wind velocity profile to trajectory' 
prediction. A less elaborate alternative is to adopt a wind profile model such as the 
Ekman spiral for wind speed approximation [Ref 5[. Since the change in the predicted 
range is not expected to be large, the wind velocity is treated as constant throughout this 
thesis. 
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II. PROPOSED C\V RADAR DESIGN 



A. RADAR SYSTEM SCHEMATICS 

The receiver design of a radar is critical to its range performance. Since no specific 
information for any C\V radar tracking system is available, the following receiver con- 
figuration is proposed. Here Figure 1 shows the schematics of the radar. 




ANTENNA 

BEAM 

INFORMATION 



Figure 1. Schematics for Proposed CW Radar Design. 
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The received signal at the frequency (/^ +X) the receiving antenna is fed to a 
narrow-band, low-noise RF pre-amplifier with a bandwidth of 1 MHz centered at 10.475 
GHz. The noise figure of such a pre-amplifier can be as low as 1.5 to 2 dB, with a 
minimum gain of 28 dB [ Ref 6 ]. This low-noise pre-amplifier can limit the receiver 
noise figure ( f„ ) to about 2 dB. The 1 MHz bandwidth is chosen according to the es- 
timated maximum doppler frequency and the fact that the target is always flying 

away from the radar. The output signal from the low-noise pre-amplifier is reduced to 
the base band at the mixer, and then passed through the video filter of bandwidth 5,. 
The filter output, at the doppler frequency/,, can be recorded either with an analog or 
a digital recorder. The bandwidth of this recorder should be wider than to allow 
truthful recording of the raw data for later analysis. The doppler signal is then fed into 
the square-law detector. 

The output of the square-law detector oscillates at a frequency of 2/. Since this 
frequency is proportional to the velocity of the projectile receding away from the radar, 
it can be integrated to obtain the projectile range information. 

The integration is most easily accomplished by counting the number of cross-over 
points of a pre-determined level by the square-law detector output. As to be explained 
later, this number, denoted by MN, is the displacement of the projectile away from the 
radar measured in half the transmitter wavelengths. 

The square-law detector, a counter, and a clock will be adequate for measuring the 

range and the velocity of the projectile receding from the radar. First, a measurement 

base, .MB, is set in terms of the number .VIN of transmitter half wavelengths. Then the 

time T„ for the output of the square-low detector to run through .VIN pre-set level 

W B 

crossings is measured. Here VIB is the radial range increase of the projectile and ■■ " •■■■■ 

■* n 

is the radial velocity of the projectile. 

As discussed in the introduction, the trajecton.' of the projectile is predicted before- 
hand. The elevation and bearing angles of the radar are updated at the end of every T„ 
interval. This information is sent to the servo control system to keep the antenna on the 
projectile. 



B. THEORY 

The theory of operation for the proposed C\V radar system will be described. The 
doppler effect is used in a C\V radar for moving target detection. This effect was found 
by Christian Doppler in 1842. It states that the frequency of a harmonic signal will shift 
when reflected by a moving target. The change in frequency is proportional to the ve- 
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locity of the target. Let denote the doppler shift, V, the radial velocity of the target 
with respect to the radar, c the velocity of light, and f, the transmitted carrier frequency. 
The doppler frequency shift is 



/a =2-^ (2.1) 



The doppler is zero when the trajectory of the projectile is perpendicular to the radar 
line of sight. It is positive for an approaching target (see Figure 2). 



Amplitude 




Amplitude 




T 

Amplitude 








^0 /q ^0 

Frequency Frequency Frequency 

(^3) {^) ic) 



Figure 2. Spectra of Received Signals: (a) no doppler shift, no relative motion (b) 
approaching target (c) receding target. 



The frequency of the signal out of the square-law detector is 2f^. Thus the number 
of crossings, MN. of this signal over a pre-set level within a time interval T„ is : 

.U.V=r,x^ (2.3) 

= r,x2-^ (2.4) 

According to equation (2.1), the radial displacement of the projectile over this period 
(Figure 3 on page 1) is : 
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Figure 3. Measurement Base and \’elocity vs. Time. 



MB = X i; 



= A/.Vx-^ 



(2.5) 

( 2 . 6 ) 



For example, if the transmitting frequency / of interest is 10525 MMz, we can set 
the number MX to 100. Then 



100 X 2.997 X 10^' 
2 X 10525 X 10^ 



(2.7) 



= 1.4234 m (2.S) 

This is a small change in target location for updating the radar pointing direction. 
Since MB is the displacement of the target along the radar pointing direction, this dis- 



7 



placement can be projected onto the ground. The ground range can be derived by add- 
ing up the individual projected ranges. 
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111. PROJECTILE TRAJECTORY PREDICTION 



To analyze the range performance of the radar, information about the trajectory of 
the projectile is required. There are many factors such as air drag, gravitational force, 
projectile shape and spin. Earth rotation,. ..etc., which can affect the trajector}'. As dis- 
cussed in the INTRODUCTION, the projectile spin mainly provides a stabilizing 
mechanism to the projectile and is assumed to have little effect on its trajector}'. The 
direction of the projectile a.xis (or heading ) is assumed to be along the direction of the 
air velocity of the projectile. Thus only the motion of the center of mass of the projectile 
needs to be considered. The projectile center of mass and the center of air pressure are 
assumed to be the same. The trajectoiy is thus reduced to that of a particle trajector>'. 

Wind velocity can increase, decrease or deviate the trajectorv' through air drag on 
the projectile. In this thesis, it is assumed that wind velocity is constant over the course 
of the projectile; it is usually called the ballistic wind. 

Drag force depends on the air velocity of the projectile. It points against the 
projectile heading. Its magnitude is proportional to the square of projectile air speed, 
the local air density, the projectile cross-section, and the drag coeOicient. Air density 
depends on the temperature profile of the atmosphere. The drag coefficient is very 
complicated and, in fact, varies with the projectile air speed. Measured data of the drag 
coefficient are available. The drag force on the in-flight projectile is given below : 

7 = -h-^ACo{V- iTy (3.1) 

where If' is wind velocity. I ' is the velocity of the projectile, Q is the drag coefficient, 
A is the cross-section area of the projectile, and p is the local air density. Here h is the 
heading of the projectile, which is a unit vector in the direction of I'— jJ'. 

The gravitational force is the only other factor affecting the trajecton.' in an inertial 
frame of reference. In the accelerating reference frame fixed on the surface of the Earth, 
there is the fictitious Coriolis force. It causes an apparent trajector}' drift. As observed 
on the ground, this drift can reach 3.1 km for a trajectory of 27.4 km ground range. 
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A. THE TRAJECTORY IN A FIXED INERTIAL FRAME 

In order lo simplify the equation of motion, the Earth-centered inertial frame (r, 
d, 0) is used. The plane containing the radar and the spin axis of the Earth when the 
projectile is fired at t = 0 is the z-x plane. The radar is moving with the Earth at an 
angular velocity of Q . Its location in this inertial frame is (>•;<. 0^, Qt) . 

The components of the projectile velocity in this reference frame are : 



7= i; 


(3.2) 


>e=v^ 


(3.3) 


>•( sin 0)<^ = J 


(3.4) 



The Lagrangian formulation of the principle of least action is convenient for deriv- 
ing the equations of motion of the projectile [Ref 7]. The Lagrangian of the projectile 
in the inertial spherical coordinate system [Ref 1 , 8] ( r, 6, (f) ), without considering the 
drag force, can be expressed as : 



L= T(r, i\r,e,4)) 



(3.5) 



L = — {)■ + r 6 + r 4> sin d) — mga( 1 — — ) 



(3.6) 



where g is the gravitational acceleration at ground level, 
equation of motion is ; 

d cL ^ sL 

di cr cr 

r = — g + r6^ +r4>" s\r\^ 6 
r 

For the 6 component : 

d cL ^ _cL_ 
di 66 

-^{r^6) = sin 0 cos 0 
di 



For the r component, the 

(3.7) 

(3.8) 

(3.9) 
(3.10) 
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and, for the 0 component: 



d cL 

di ^ 



8L 

c4> 



-|-(-^(^sin'0) = O 



(3.11) 



(3.12) 



The nonlinear drag force can be included. The complete equations of motion become : 
d 



di 



Vr = -iTfs + 



I I I 

* P ^ * 6 



+ /r 



= cot 01'^ + ,/, 



(3.13) 



(3.14) 



— (r sin 01'^) = /-sin 04 



r = a + /2 



(3.15) 



(3.16) 



where h is the height from the Earth surface and a is the radius of the Earth. The ratio 
of projectile height to the Earth radius can be expressed as ; 



-L = 1+A 

a ^ ' a 



(3.17) 



The projectile velocity is : 

V = v;+r^e+rj 



(3.1S) 



= >r + >-00 + r sin 



and the wind velocity is : 



(3.19) 



>r = jVgd + + u\r 



(3.20) 



Note that are the components of the drag force divided by the mass of the 

projectile : 



7 = 



(3.21) 
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(3.22) 



= Lr +fee +fj 



The equations of motion can be discretized. Assume that the gravitational and the 
drag force acceleration are approximately constant over the period T„ when the projectile 
moves from r, to the equations of motion are approximately : 



'n+] 



T„ 



+ 






(3.23) 



- r;e 



(3.24) 



f = 

* n 

From the above equations, if r„. 6^. </>„ and r„. </>„ are known, then r„.j, 

and can be derived as follows : 

With measured 7),. equation (1) provides , then 



(3.25) 



= fr. + 2 

Equations (3.26) and (3.24) give then 



• n+l 



= + 






T 

* n 



Equations (3.27) and (3.23) give then 



(3.27) 



4>n+\ = + 



+ <^n+\ 



T 

* n 



(3.2S) 



Note that the muzzle location is and the muzzle velocity of the projectile, which 
is assumed to be known, is ! Through the above procedures, the full trajectory of the 
projectile can be constructed by computer simulation. 
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B. GROUND-OBSERVED TRAJECTORY 

In order to obtain the trajector}' as observed from the ground, the position of the 
projectile should be converted into its relative location from the radar antenna and from 
the gun. The necessarv' coordinate transformation matrices are given below (Figure 4 
on page 13). 




Figure 4. Inertial Coordinate Systems. 

1. Trajectory Relative to Radar 

To set the elevation and bearing angles for the radar and to estimate the doppler 
shift, the location and velocity of the projectile relative to the radar as viewed from the 
ground are needed. On the ground, the projectile motion is most easily described in an 
inertial frame in which, at time, t, the radar is passing through the z' —V plane instan- 
taneously. The z' axis is chosen so that it coincides with the z axis. The transformation 
from the fixed inertial spherical coordinate system to this (.v',y' ,z') system can be ac- 
complished through the following rotation matrices , and R^^, : 
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R,{0.i) = 



cos 12/ sin 0.1 0 



(3.29) 



— sin Oi cos Oi 0 
0 0 1 



Rreci<i>p) 



COS 4>p 0 — sin 4>p 

sin (^p 0 cos 4>p 

0 1 0 



Rcy^Op) = 



sin 6p 
cos dp 
0 



cos Bp 0 
— sin Bp 0 
0 1 



(3.30) 



(3.31) 



a. J'elocity Transformation 

The projectile velocity in the (.v', y'. z') coordinate system can be obtained 
by first changing the results in section lll.A from the spherical coordinate system into 
the rectangular coordinate system (,\, y. z) ; 









^p.v 


= Rr,fd>p)Rcyl{Bp) 








^ pc 



Since the (.v', y', ;') coordinate system is obtained by rotating the (.\, y, z) 
coordinate system by an angle Oi around the z-a\is. 







’’p/ 


'PT 


= /?,(f2/l 


^ py 









(3.33) 



= Rreci<t^p-0,)R,p^Bp) 



r. 



p'- 



pe 



pc 



The velocity of the radar, I r, in the {x', y' ,z') coordinate is : 



(3.34) 
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(3.35) 



^ /ex' 




0 


J'/ev' 


~ ■^m(0)^c>'/(^/e) 


0 









sin 0/j 


cos 


o' 


0 


0 


0 


1 


0 


cos 


— sin 


0 


^ R 4 > 



0 

= 

0 _ 

Hence 

b. Projectile Location Relative to the Radar 

The relationship between {x'.y', z') and (r, 6, 4>) is given by ; 
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(3.36) 



(3.37) 



(3.3S) 



(3.39) 



(3.40) 



(3.41) 



(3.42) 
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c. Radar Coordinate 



Substitute {r;^, 0^,. Qt) into equation (3.42), 
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2. Trajecton Relative to Gun 

Relative to the gun (muzzle), a similar inertial frame (.v " ,y “ , z " ) in which the 
gun is instantaneously passing through the z " - .r " plane at time t can be used. The 
z “ a.xis coincides with the z-axis. The coordinate transformation is : 



X " 




/• sin 6 cos(^ — 4>c ~ 




= 


r sin 0 sin(<5 — do — ^0 






r cos 6 



(3.44) 



The firing plane of the gun is the plane containing its barrel and passing through 
the center of the Earth, for artillery, it is customary to use the coordinate system 
(.v,._iy z^) in which the gun muzzle is the origin and the z,— Xf plane is the firing plane. 
Here z, is the vertical direction and x, is along the horizontal from the gun. The relation 
between (.w.yv. z.) and (.v " . " . z ") depends on the gun bearing angle Q which is the 

angle of the .t, axis measured counterclockwise from the East. If the the gun is oriented 
along is located at I'c . (f)^. the vector — fe in the coordinate (.xv.yv. z,) is : 
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Note that /. - is the relative location of the projectile from the gun. Hence 
its components can replace (.r " , y " ,z " ) directly in equation (3.45). 
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IV. RANGE ANALYSIS AND POWER REQUIREMENT 
A. RANGE EQUATION 

The power requirement of the CW radar can be obtained from the radar equation 
[Ref 9] : 



P, = 






c,a):r 



(4.1) 



In equation (4.1). the radar cross-section {o). the maximum range (R), and the 
propagation factor (F) will be discussed in this chapter in detail. The Boltzman constant 
k is 287.05 n-m kg-°A'. The absolute temperature T, is 290°A' . The wavelength / is 

0.02S5m. Other parameters are based on the configuration of the CW radar as proposed 
in Chapter 2. They are defined as follows. 

1. The Receiver Noise Figure 

The receiver noise figure (AJ is the ratio of the signal-to-noise power ratio at the 
receiver input to the signnl-to-noise power ratio at the receiver output, assuming that the 
receiver temperature is 290' A' . The A„ value for the receiver considered is assumed to 
be 2 dB. 

2. System Loss Factor 

The loss factor consists of the fractional loss of transmitter power delivered to 
the transmitting antenna (L,). the fractional loss of the receiving antenna output power 
delivered to the receiver input (A,), and miscellaneous loss (A„). They are discussed be- 
low. 

a. Transmission Loss 

The transmission loss is the loss in the transmitting or the receiving portion 
of the transmission line. Since the transmitter is light-weight and is to be mounted very 
close to the antenna. L, is estimated to be about 1.2 dB due to the use of narrow band- 
width rotarv' joints. The receiver pre-amplifier can be mounted directly to the receiving 
antenna output, followed with the mixer immediately. This arrangement will reduce A. 
to about the same level as A,. 
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b. Miscellaneous Loss 

The miscellaneous loss (L„) for this CW radar system is mainly due to the 
quantization loss at the comparator. Since only two-level (crossing) comparison is 
made, it is not expected to be significant and will be ignored. 

Base on the above argument, the total loss is estimated at 2.4 dB. 

3. Signal-to-Noise Ratio 

For a tracking radar, it is important that the receiver output signal is sufficiently 
high above the noise level so that information carried by the signal can be extracted re- 
liably. With the proposed receiver configuration, the signal-to-noise ratio (S N) of the 
output from the video filter is the figure of significance. First of all, this signal is re- 
corded for later analysis. Its minimum S N has to be assured. In real-time, after passing 
through the square-law detector, this signal generates the timing sequence for updating 
the radar pointing direction and for evaluating the projectile radial velocity. It drives the 
gain control to keep the reference signal at a preset level. It is also the input into the 
comparator following which the number of signal cross-overs is tallied and the elapsed 
time is clocked as explained in section ll.A. 

The square-law detector and the comparator perform threshold detection. The 
reference signal acts as the threshold. 1 he comparator outputs a + 1 signal when the 
detector output is above the reference level; it outputs a -1 signal when the detector 
output falls below the reference. The reference signal level, which is set at a fraction of 
the detector output signal, determines the accuracy of the real-time output of the system. 
In what follows, the minimum S N ratio and the threshold setting in terms of the 
threshold-to-signal power ratio are discussed. 

The signal input to the comparator is oscillating at twice the doppler frequency. 
Within each cycle of the signal, the portion rising above the threshold level to the peak 
and back to the threshold level can be considered as a pulse. This pulse is contaminated 
with noise, which is assumed to be white. The probability of each contaminated pulse 
falling under a certain threshold level is well-known. This probability of missed de- 
tection is the fraction of signal pulses failed to be detected by the comparator. Between 
two pulses, the number of false alarms is the number of extra cycles counted. These two 
effects tend to cancel each other so the optimum threshold setting is achieved by equal- 
izing the fraction of missed pulse and the fraction of noise pulses. The program to 
compute these fractions are listed in Appendix D (Ref 10]. 

Figure 5 on page 19, Figure 6 on page 20, Figure 7 on page 20. Figure 8 on 
page 21, and Figure 9 on page 21 show the fractional errors due to missed pulses and 
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noise pulses respectively. If the signal-to-noise ratio is 10 dB, the comparator output 
error is 5.55% when the threshold-to-signal power ratio is set at 0.455. Stronger signal 
allows the threshold to be set at a lower fractional value with an improved error ratio. 
For example, for a signal 12 dB above the noise, the threshold setting of 0.38. or -4.2 
dB below the signal, will limit the error to I.151o. In this thesis, the minimum signal- 
to-noise ratio is chosen to be 10 dB. 
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Figure 5, Missed Detection and Noise Error vs. Threshold-to-Signal Po»er Ratio 
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Figure 6. Missed Detection and Noise Error vs. Threshold to Signal Po^^er Ratio 
S N= 10.5 dB. 




Figure 7. Missed Detection and Noise Error vs. Threshold to Signal Po«er Ratio 
S N= 11 dB. 
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Figure 8. ^^issed Detection and Noise Error vs. Threshold-to-Signal Po>\er Ratio 
S N= 11.5 dB. 
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4. Transmitting and Receiving Antenna Gain 

Separate but identical transmitting and receiving antennas are used to reduce 
direct coupling of the transmitter power into the receiver. Since this separation is small, 
it is neglected in the analysis. In x-band. antennas with one-degree beamwidth are easily 
achievable. Based on this beamwidth. the antenna gain (G) [Ref. 9] is : 



G = 



27000 

AO- 



(4.2) 



= AA.MiB 


(4.3) 


and the diameter D of such an antenna is 






(4.4) 


= 72 X " 


(4.5) 


= 2 m 


(4.6) 



5. Video Bandwidth and Filters 

The required minimum video bandwidth {BJ is a function of the doppler fre- 
quency. .-\s previous mentioned, the bandwidth should be greater than . Although 
the doppler frequency/, of the projectile increases rapidly right after firing, it decreases 
monotonicaliy aftertvards by a factor of more than 4 over its trajecton.’. An example is 
shown in Figure 10 on page 23. 




Referring to the range equation (4.1), in order to use the lowest transmitted 
power possible for full range tracking, the radar range and video bandwidth can be 
traded off Note that, assuming a constant radar cross-section, the transmitter power 
depends on the product A filter of narrower bandwidth may be employed when 

the projectile reaches a greater distance. Adoption of a set of three switchable doppler 
filters is recommended. The following is an e.xample for its implementation based on the 
doppler frequency of the projectile given in Figure 10. 

The ma.ximum doppler shift from the projectile shown in Figure 10 is a 
little less than 54 kHz. The maximum range R is about 27 km. Three filters are to be 
used over three different ranges. Each of the filters will have a bandwidth of 1.2 times 
the maximum expected doppler frequency over the range in which it is to be used to al- 
low for a 20°o estimation error. The transmitter power level will be minimized if the 
fourth power of the maximum range times the maximum doppler shift in the corre- 
sponding range is a constant. 
a. First Filter 

When the projectile is within 16.06 km slant range, the video bandwidth is 
chosen to be 1.2 x Thus : 
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B,= 1.2 



(4.7) 



= 63.9 A//r 



(4.8) 



b. Second Filter 

When the projectile moves beyond 16.06 km. but within 23.73 km slant 
range, the doppler frequency has a maximum value of about 22.03 kHz. The video 
bandwidth is : 



B, = 1.2 X 22.03 kHz 



(4.9) 



= 26.44 A//c 



(4.10) 



c. Third Filter 

When the projectile reaches 23.73 km slant range, the maximum doppler 
frequency is about 12.7 kllz. The video bandwidth is : 



B, = 1.2 X \1.1 kHz 



(4.11) 



= \ 5. 2 kHz 



(4.12) 



Since is now less than 1 4 of the original value when only one video filter is used, 

the minimum required transmitter power will be reduced by more than a factor of 4. 
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B. TRAJECTORY SIMULATION AND RESULTS 



Based on the procedures for trajecior>’ estimation presented in Chapter 111. a simu- 
lation program was written and is included as Appendix A. It was written in PC-Matlab 
3S6. The program is described as follows. 

I. Input and Parameters. 

a. Earth 

The operation is assumed to take place in Taiwan. The radar is located at 
the longitude of 121'’36'27 " East and the latitude of 24°36'3S " North. The Earth ra- 
dius is approximately 6.3677225 x 10*^ m (average value of the artic radius 
6.356912 X 10'' m and the equator radius 6.37S533 x 10* m). The gravitational acceler- 
ation at sea level is 9. SI m/ sec% The angular velocity of Earth spin is 
7.2921 1S44 X 10--‘ rad sec (Ref 1 1]. 

b. Radar 

The transmitting frequency is 10525 MIlz. The measurement base MB can 
be varied to change the update rate of the antenna pointing direction. For both the re- 
ceiving and the transmitting antenna, the beamwidth is T and the gain is 44.3 dB. The 
antenna height and the distance between the gun and the radar can affect the initial 
readings of the doppler signal. 

c. Gun 

The firing elevation angle is assumed to be 45° and the azimuth is counter- 
clockwise from the East and is also assumed to be 45°. Both the elevation angle and the 
azimuth angle can be varied. The distance between the gun and the radar can be ad- 
justed. (Variations of up to 2000 m are possible at the test site.) The barrel length of 
the gun under consideration is 7 m. 

d. Projectile 

The target is a 155 mm calibre artillerv' projectile. Its weight is 45.4 kg. with 
a total length of 93S mm. The maximum muzzle velocity is 850 m sec. 

e. Environment 

n j Temperature and air density profdcs. The temperature at sea sur- 
face is assumed to be 288° A' and the temperature at the projectile height h is based on 
the following profile [Ref 12] : 

Th = (4.13) 
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where h is the height from sea level, Oj- is the temperature lapse rate (-0.0065 'K m). 
The temperature (T,) is used for calculating the air density (pj at the projectile: 



^ O 



(4.14) 



In equation (4.14) p, is the air density at sea level (1.225 kg nT). and R is the gas con- 
stant (2S7.05 n-m kg-°A'). 

Wind profile. The wind is assumed to be constant during the flight 
time over the projectile trajectory. The value used is the so-called ballistic wind. It can 
be expressed as : 

W = Wf + WgO + W^i) (4.15) 



In this thesis, the wind is 6 m sec, North-West. It should be noted that the extent to 
which the estimated trajectory’ can be modified by the actual wind profile is subject to 
further study. 

Mach number. 1 he Mach number is the ratio of the projectile air 
speed to local sound speed. The sound speed at the height of the projectile is : 

a, = (4.16) 

where K is 1.4 for idea gas. The Mach number (.M ) is given by : 

I f' - » k I 

M = (4.18) 

where | (' — H' | is the projectile speed with respect to air. 
f. Drag Coefficient 

The drag coefficient {Cf) is required for calculating the drag force. It de- 
pends on the shape of the projectile. The Q curve is measured experimentally. It is 
given as a function of Mach number M. For a 155mm projectile, the Cp curve is shown 
in Figure 11 on page 27. 

2. Results of Simulation 

The results of trajectory simulation are described below. The variable names 
used in the program are given in the parentheses when they are first mentioned. 
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a. Radar Observed Trajectory 

The slant range (SRR) is the distance from the radar antenna to the in-flight 
projectile. The curve showing the projectile height (Height) versus its slant range from 
the radar is shown in Figure 12 on page 2S. The total flight time is 87.4 seconds, and 
the maximum range is 27.4 km. 

b. Gun Based Trajectory 

The ground range is measured from the gun to the projectile. The projectile 
height vs. in-flight time (t) and the projectile height vs. ground range are shown in 
Figure 13 on page 28 and Figure 14 on page 29. The maximum ground range (Grange) 
is 27.4 km. 




Figure 11. Drag Coefficient vs. Mach Number for 155 mm Projectile. 
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Figure 12. Projectile Trajectory: Height vs. Slant Range. 




Figure 13. Projectile Trajectory: Height vs. Time. 
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The cross range (CRange) measures the drift of the projectile from the firing 
plane. The firing plane is the vertical plane containing the gun barrel. In the case con- 
sidered, the maximum cross range is about 6.3 km South-East of the firing plane. The 
point of impact is about 3.14 km away from the firing plane. The curve of cross range 
vs. ground range is shown in Figure 15 on page 30. 
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Figure 15. Projectile Trajectory: Cross Range vs. Ground Range. 
c. Radar Tracking Parameters 

The azimuth (traz) and elevation (trel) antenna pointing angles are required 
for tracking. The curves of these angles vs. in-flight time when the gun-radar distance 
is 500 m are shown in Figure 16 on page 31 and Figure 17 on page 31. 

The azimuth and elevation angular speeds of the antenna are shown in 
Figure IS on page 32, and Figure 19 on page 32. 
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Figure 16. Radar Azimuthal Tracking Angle vs. Time. 




Figure 17. Radar Elevation Angle vs. Time. 
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Figure 18. Radar Azimuthal Tracking Angular Speed vs. Time. 




In-flight time (sec) 



Figure 19. Radar Elevation Tracking Angular Speed vs. Time. 
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The antenna servo system has to point the antenna to within half a beamwidth in the 
direction of the projectile. Thus the half beamwidth divided by the angular speed of the 
antenna is the required response time of the servo system. The servo response time vs. 
the projectile in-flight time is shown in Figure 20 on page 33. From the above results, 
it can be seen that over the first few seconds, the elevation angular speed is ver\' high. 
The distance between the gun and the radar can be increased to reduce this speed. 




Figure 20. Servo System Response Time vs. Time. 



d. Doppler Frequency and Velocity 

The doppler frequency (/Q) is a directly measurable data by a CW radar. 
The doppler cycle can be detected with the square law detector. In the proposed CW 
radar configuration of this thesis, the doppler cycle is derived from the time interval 7'„ 
for each measurement base traveled. This timing sequence has to be simulated in the 
tracking program. In the simulation program, the doppler velocity (V^^^,) is calculated 
from the projectile velocity relative to the radar and the radar pointing direction. Then 
T„ is estimated from the time required to travel the distance MB with this doppler ve- 
locity. The curves of doppler frequency vs. slant range, doppler frequency vs. time, and 
doppler velocity vs. time are shown in Figure 10 on page 23, Figure 22 on page 34, and 
Figure 21 on page 34. 
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Based on the above results from trajectory' simulation, the projectile radar 
cross section and the transmitter power can be estimated. They will be discussed later. 



xio^ 




Figure 21. Doppler Frequency vs. Time. 




Figure 22. Doppler \'elocity vs. Time. 
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C. PROJECTILE RADAR CROSS-SECTION 



The projectile radar cross-section (RCS) is denoted by a. It is the area intercepting 
the amount of plane wave power which, when scattered equally in all directions, 
produces an echo at the radar. In other words, the RCS is given by : 

o = 4rtr\-^f (4. IS) 



where is the scattered electric field strength from the target, r is the distance from the 
target in its far field where the quantity £, is measured, and £, is the incident plane wave 
field strength at the target. 

Determination of the RCS of a complex target is difficult. Usually the radar cross- 
section of a target depends on the following factors: a) target size, b) target shape, c) 
radar frequency, d) target aspect angle, and ej radar polarization. The RCS is greatly 
simplified if the target size is large compared to the radar wavelength. Consider the ra- 
dar cross-section of a simple sphere. The RCS as a function of the circumference of the 
sphere measured in wavelengths {2nal'/.) is shown in Figure 23 on page 36. The region 
where the target size is large compared to the radar wavelength, i.e.. where 2nal/. is 
greater than 10. is called the optical region [Ref 9). The RCS in this region is roughly 
the physical cross-section at area of the target. 

The 135 mm projectile is in the optical region at the transmitter frequency 
(10525 MHz) of interest. Its RCS value will be approximated with its geometrical 
cross-sectional area. Because of the rotational symmetry of the projectile, the radar 
cross-section of the projectile is a function of the angle between the radar pointing 
direction {£) and the projectile heading (/i). 

The radar pointing direction D is given by : 



V P — y p 

D = -4- — r 

’ - >R 



(-1.19) 



where 7p is the projectile position vector, 7p is the radar position vector as defined in 
Chapter III. 
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The heading of the projectile, h, is in the direction of its air velocity. It is also de- 
fined in Chapterlll: 



In terms of h and D, the aspect angle of the projectile as seen from the radar is 
e,,^ = cos-'(D./i) (4.21) 



Note that 0 < < n. For simplicity, the projectile is considered as a circular cone on 

top of a circular cylinder, as shown in Figure 24 on page 37. The height of the cylinder 
is 414.4 mm. The height of the cone is 523.6 mm. The diameter of the cylinder and the 
base of the cone is 155 mm. 
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The vertex angle of the cone, 2a, is given by; 



a = tan 



-1 155 



2 X 523.6 



= 0.1469 rad. 
= 8.42° 



(4.22) 



(4.23) 

(4.24) 



The physical cross-section of the cone is the projected area of the projectile along D. the 
radar pointing direction, onto a plane perpendicular to D. This cross-sectional area will 
be discussed below according to whether | Z) • /i | < cos ex or | Z) • /) | > cos a. 

1. Case 1 

As \ D • h \ > cos a, the projection of the projectile along D is shown in 
Figure 25 on page 38. 
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The loial RCS value is 

c,= rM-\b^h\^lab\b^h\ (4.25) 

where a in the radius of the cylinder and b is the height of the cylinder. 

2. Case 2 

As \ b>h \ <co$y., the total RCS is given hy equation (4.20) (Figure 26 on 
page 39). 
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Figure 26. Radar Cross Section Case 2. 



<7; = <7] + CT2 + 03 + 04 



(4.26) 
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(4.27) 



(4.28) 



(4.29) 



03 = \ b X h \ 



(4.30) 
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(4.31) 



, \b-h\ 

Oi = tm' r 



r ^ r\ 

A D cos‘a — h cos 

A’ = — — ■ ^ =- (4.32) 

sin cos a + sin a cos 

where L is the height of the cone and .V is the normal vector of the surface A,. The de- 
tailed derivation of the RCS equations is given in Appendix B. 

From the above results, the RCS value of the projectile in-flight can be pre- 
dicted. It is shown in Figure 27 and Figure 28 on page 41. From above figures, it can 
be seen that the maximum RCS is about 0.09 m- and the minimum RCS is about 
0.02 /7I-. Since the aspect angle of the projectile changes ver>' fast right after it is fired, 
the RCS value varies rapidly at the beginning of the trajectory. 




Figure 27. The In-Fligl t Projectile Radar Cross Section vs. Aspect Angle. 
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Figure 28. In-Flight Projectile Radar Cross Section vs. Slant Range. 
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D. PROPAGATION EFFECTS 



Propagation of radio waves through the Earth's atmosphere is significantly affected 
by the propagation medium. The presence of the Earth complicates the situation in 
three ways: (a) the surface of Earth scatters radiation and produces an interference pat- 
tern, (b) the Earth casts a shadow and gives rise to diflraction phenomena, and (c) the 
Earth's atmosphere is inhomogeneous and produces refraction effects [Ref 13]. It is 
customary to group all these factors into a single quantity called the propagation factor 
(E). In this thesis, the Engineer's Refractive Effects Prediction system (EREPS) [Ref 
14] is utilized to compute the propagation factor based on local environment. 

EREPS uses ray tracing for radio wave propagation within the horizon. It includes 
the bending effect of the refractive inde.v variation of air and the divergence factor when 
the ray is reflected from the Earth. The field at the target is the sum of the direct ray 
and the reflected ray. resulting in interference peaks and nulls which become significant 
when the target falls near the horizon (Figure 29). 



^target 




Figure 29. Reflection and Interference from a Spherical Earth. 
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For the over-the-horizon region, the waveguide mode approach is adopted by 
EREPS. The Earth's atmosphere is modeled as a leaky waveguide. The Earth's surface 
is the lower boundary of the guide. .VI odes of high attenuation along the horizontal di- 
rection are called leaky modes. .Modes of low attenuation are called trapped modes. In 
x-band. the effect of the evaporation duct is significant. Trapped modes occur when the 
evaporation duct height is large or at least comparable to the radar wavelength (see 
Figure 30 on page 44) (Ref 15). 

1. Computation 

A typical EREPS output is given in Figure 30 on page 44. The local annual 
average meteorological conditions for propagation are available [Ref 13). The following 



parameters are used : 

1. Evaporation duct height : 13.1m 

2. k (effective Earth radius factor) : 1.45 

3. Wind speed : 6ms 

4. Absolute humidity : lOg m^ 

5. Surface refractivity : 370. 



Note that the projectile is at a particular height when it arrives at a range. This 
computation has to be repeated for each point along a projectile trajectory. This is done 
for several trajectories at various ranges (Figure 32 on page 45 and .Appendix C ). Also 
note that, in Figure 32 on page 45. the ground range is given a negative sign because it 
is measured from the point of impact toward the radar. 

2. Summary of Results 

a. Signal Strength Enhancement 

From Figure 32 on page 45. it can be seen that, with ducting, which is 
present most of the time at the test site, the return signal will be substantially enhanced 
from the corresponding free space value. (Note that the computation is carried out for 
a target at 1 m height, i.e., when the projectile is near its point of impact. The applicable 
range for the projectile of interest lies between 27 to 28 km.) Thus ducting is helpful in 
enhancing the quality of the recorded full range tracking data. 

b. Strong Interference 

As shown in Figure 32 on page 45, the return signal will exhibit strong in- 
terference when the projectile reaches about 250 m within the point of impact. The more 
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than ± 15 dB fluctuation makes it diflicult to track the projectile in this region using the 
simple CW system proposed in this thesis. On the other hand, since this strongly fluc- 
tuating region is within a few hundred meters from the impact point, this small portion 
of the trajectory can be estimated with good accuracy. Hence it is desirable to select a 
'stop-processing point' at which the CW real-time radial velocity and range computa- 
tions will be halted. The radar antenna will be guided by the predicted trajecton.', and 
the return video signal will still be recorded. 
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Figure 31. Propagation Factor nitli 13.1ni Duct Height. 




Figure 32. Propagation Factors vs Ground Range from Impact Point. 
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3. Stop-Processing Point 

It should be noted that EREPS does not take into account the antenna radi- 
ation pattern. Since the antenna to be used have a beamwidth of 1° only, interference 
between the direct ray and the reflected ray will not be important until the antenna ele- 
vation angle is lowered to within 0.5° above the horizontal direction (Figure 33 on page 
46). 




This criterion can be used to determine the stop-processing point as follows; when the 
system records an antenna elevation angle falling below 0.5°, the square-law detector 
and the comparator will be switched off. The servo will steer the antenna according to 
the predicted trajectory smoothly toward the projectile point of impact. Characteristics 
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for two such stop-processing points for a radar-to-gun distance of 500 m and two dif- 
ferent radar heights of 10 m and 30 m respectively are displayed in Table 1 on page 47. 
Note that these point are close to where the propagation factor starts to fluctuate 
strongly as predicted by EREPS. 

Since the processing of real-time data is stopped before contributions from the 
reflected ray become important; in what follows, the propagation factor will be consid- 
ered as unity. 



Table 1. STOP PROCESSING POINT. 



Antenna height (m) 


10 


30 


Radar to gun distance(m) 


500 


500 


Projectile height (m) 


308 


305 


Slant range (km) 


27.53 


27.53 


Impact range (km) 


27.69 


27.57 


Ground range to impact point (ni) 


15S 


157 


Elevation angle (degree) 


0.49 


0.45 


Elevation angle at impact point (degree) 


-0.15 


0.19 


Elevation angle diflerence (degree) 


-0.64 


-0.64 


.Azimuth angle (degree) 


36.5 


36.5 


.Azimuth angle at impact point( degree) 


36.9 


36.9 


.Azimuth angle difl'erencc to impact point(degrce) 


0.39 


0.39 


In-flight time (sec) 


S6.4S 


S6.4S 
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E. POWER REQUIREMENT 

From equation (4.1). the power requirement for tracking up to the stop-processing 
point is shown in Figure 34. In order to reduce the required transmitter power re- 
quirement, a multiple filter design is proposed for this system. 

The multiple filter design is based on the characteristics that, after an initial rise, the 
doppler frequency decreases with increasing slant range (Figure 35 on page 49). By 
switching to video filters of narrower bandwidths at longer ranges, the coherent inte- 
gration time of the radar is increased hence the power requirement can be lowered. 
Combining the fixed range independent parameters in the radar range equation, the 
power requirement can be expressed as : 



B 

P, = (3.7339 X 10"*") X ( — ^ )n^a^ 



(4.33) 



= C X ( 4- 



(4.34) 



C = 4.4S07 X 10“" 



(4.35) 



where 5. is \.2f^. 




Slant range ( n ) 



Figure 34. PoMer Requirement for Full Range Tracking nithout Multiple Filters. 
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Since f^. a and R are functions of the trajector>' (Figure 35 on page 49, and Figure 36 
on page 50), if the tracking range is divided into three regions 1, 2, and 3, and one filter 
is used for each region, the transmitter power requirement becomes : 



R 



p, = C(-^ (F,) 



1 .max 



~ a h.max (fdh,max 



R 



— a h,max ifdh.n 



where the constant C is the same as that of equation (4.34). 



xiCH 




slam range (m) ^10" 



Figure 35. Doppler Frequency vs. Slant Range. 



(4.36) 

(4.37) 

(4.38) 



49 




Figure 36. Radar Cross Section vs. Slant Range. 

The power requirement written in terms of that of the single filter case is; 

(4.3 

where C^. the power reduction factor, is the factor by which the power requirement is 
reduced. 

In order to find the optimal switch-over points for the filters, the logarithm of f, 
versus R^lo, normalized to and respectively, is plotted (Figure 37 on page 

51). 

.Assume that the first switch-over point is located at (-v,.yi). and the second one is 
located at (.r^,.!-). These values can be determined from Figure 37 on page 51 by ob- 
serving that, at the optimal switch points : 

(4-40) 

= .v; (4.41) 






a f 

a )mc4xJ Jr\cix 

c 
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The optimal locations so obtained are listed in Table 2 on page 52. The power re- 
duction factor is = 4.1S6, as seen from Figure 37 on page 51. 

The power requirement as a function of slant range when three filters are used is 
plotted in Figure 38 on page 52. The antenna height is 10 m and the radar-to-gun dis- 
tance is 500 m. The required power is 37 Watts. A conventional C\V tracking system 
will track to the summit point only. This point is located at a range of region 1 and re- 
quires a transmitter power of 19.5 Watts. Almost full-range tracking is possible by 
doubling the conventional power through the use of multiple filters. 
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Table 2. MULTIPLE FILTERS SWITCH POINTS 



Parameters 


region 1 


region 2 


region 3 


stop track- 
ing point 


Slant range (km) 


16.06 


23.73 


27.53 


27.53 


Maximum X (kl U) 


53.25 


22.03 


12.72 


... 


N'ideo bandwidth 
(kHz) 


63.90 


26.44 


15.27 


... 


Radar cross section 
(m-) 


0.0385 


0.0438 


0.0862 


0.0903 


Required power 
(^^'atts) 


36.29 


36.23 


36.29 


... 
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V. CONCLUSIONS 



It is demonstrated in this thesis that full-range projectile tracking, with real-time 
track processing stopped at a range within 160 m of the projectile point of impact, can 
be achieved if the transmitter output power is 37 watts or higher and a set of three 
switchable filters are used. If a single filter is used, an output power of 160 watts is re- 
quired. Either configuration can be implemented with a contmercially available, small 
size, air-cooled RF amplifier which can be mounted near the transmitting antenna. 

The tracking relies on the predicated trajectoiy. Although real-time doppler infor- 
mation is available from the radar, it is used only to confirm that the projectile is on 
track. A algorithm to utilize this real-time information for updating the predicted tra- 
jectory' is desired. Kalman filtering appears to be a promising technique which deserves 
further investigation. 
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APPENDIX A. PROJECTILE TRAJECTORY SIMULATION AND 

RANGE ANALYSIS 



clear 
t2=clock; 
format long e; 

% The first capital letter designations are: P for the projectile; 

% G for the gun (muzzle); R for the radar, 

% The numerals 1, 2, 3 are for the radial, theta and phi components in 
% an inertial spherical coordinate system; 4 is for the radial component 
% in the cylindrical coordinate system associated with the spherical 
% coordinate. The origin is at the center of the Earth; the polar (z) 

/o axis is the axis of spin of the Earth. 

% East and South designations are local quantities in a frame rotating 
% with the Earth 



% constant declaration 
EarthT=86l64. 06; 

EAV=2'’'pi/EarthT; 
a=(6356912+6378533)/2; 
g=9.81; 

RLATI=24+36/60+38/3600; 

RL0NGI=121+52/60+27/3600; 
f=10525E6; 
lambda=3E8/f ; 

?oinput parameters 
! del W’W 
diary ww 

RB=input( ’ radar to gun base distance(m) =? 

% RBaz=input( ' radar to gun direction from east=? *)*pi/180 
RBaz=0; 

?o MN=input( ’ Measurement base number setting =? ') 

HX=1000; 

Rh=input( ’ Radar height(m)=? ') 

?o Gel=input( ’ firing elevation angle( degree)=? ^)/180*pi 
Gel=45/ 180*pi; 

% Gaz=input( ^ firing azimuth angle( degree)=? *)/180’’*pi 

Gaz=45/ 180*pi; 
diary 

answer=input( ' If you want reenter parameter please enter 1 else enter'); 
if answer=l 
return 
XX a 

end 

Bl=5; % gun barrel length 



% sidereal period of Earth spin 
% Earth angular velocity 
% approximate Earth radius 
% gravitational acceleration 

% latitude 24o36'38" (radar site) 

?o longitude 121o52’27" (radar site) 
®o radar transmitted frequency Hz 
% radar transmitted wavelength m 



A=pi*(0. 155/2)**2; 
m=45. 4; 

FAC=0. 5*A/m; 



% projectile cross section area 
% projectile mass (kg) 

% coefficient for computing drag 
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MB=MN*( lambda/2); 

% Radar Location at t=0 
Rl=a-fRh; 

R2=(90-RLATI)*pi/180; 

R3=0; 

R2s=sin(R2); 

R2c=cos(R2); 

R4=Rl^^R2s; 

Rz=Rl''R2c; 

RV=R4^^EAV; 

% Firing Direction and Speed 
PVo=850; 

Gels=sin(Gel); 

Gelc=cos( Gel) ; 

Gazs=sin(Gaz); 

Gazc=cos(Gaz); 

% Gun (muzzle) Location at t=0 
Gh=Bl*-'Gels; 

Gl=a+Gh; 

Gl=Bl“’'Gelc; 

RGs=’RB’''sin(RBaz) -Gl*Gazs; 
RGe=RB*cos( RBaz)+Gl*’'Gazc; 
G2=R2+RGs/Gl; 

G4=Gl*sin(G2); 

G3=R3+RGe/G4; 

Gz=Gl*cos(G2); 



% measurement base 
% the inertial z-x plane 

% z-coordinate of radar 
% radar velocity 

% muzzle speed of projectile 

% muzzle height 

% distance, south of radar 
% distance, east of radar 



% Displacement Transformation 
G2c=cos(G2); 

G2s=sin(G2); 

DTG=(0,GAZC,GAZS; 0, -GAZS,GAZC; 1 , 0 , 0-*-G2C , 0 , G2S; 0 , 1 , 0; -G2S,0,G2C); 



% Projectile Initial Location 
P1=G1; 

Pln=Gl; 

P2n=G2; 

P3n=G3; 

P4n=G4; 

% Projectile Initial Velocity 
PVln=PVo*Gels; 

PVg=PVo*Gelc; 

PVe=PVg'’'Gazc; 

PAV2n= - ( PVg-'-'Gazs ) / G 1 ; 
PAV3n=PVe/G4-fEAV; 



% velocity along ground 
% velocity, east 
% angular velocity, theta 
% angular velocity, phi 



% Drag Force Considerations 

% Wind Profile: assumed constant in this program. 

% The geostrophic wind should be measured and the Ekman spiral 
% model can be adopted for a better approximation if the true 
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% wind profile is not available. 

% Wind=input( * Enter the wind speed m/sec at muzzle height*) 

Wind=6; % wind velovity at muzzle height 

% WindBearing=input( * Enter the true wind bearing in degrees*) 
WindBearing=135; % clockwise from north 

Waz=( ( 90-WindBearing)+180)*pi/180; 

We=Wind-cos(Waz); % wind velocity, east 

Ws=-Wind*sin( Waz) ; % wind velocity, south 

% Temperature Profile (assumed linear) 

% Th=input( ' Enter the height (m) at which temperature is measured*) 
Th=10; 

% TREF=INPUT( ’ENTER THE MESURED TEMPERATURE IN CELSIUS*) 

Tref=15; 



TrefK=273. 16+Tref; 

Thl=Th+a; 

S =-0.0065; % TEMPERATURE LAPSE RATE K/M 

% TP=TREFK+S*(P1-TH1); % AIR TEMPERATURE AT THE PROJECTILE 



O' 

/o 

o/ 

/o 



o/ 

/o 



O/ 

/o 



AIR DENSITY 

DH=INPUT( 'ENTER THE HEIGHT (M) AT WHICH AIR DENSITY IS MEASURED') 



DH=GH; 

DA0=INPUT( 'ENTER THE MESURED AIR DENSITY IN KG/CUBIC METER') 
Dao=l. 225; % air density at sea level 
TDao=TrefK+s*(Dh-Th); 



ss=287. 05; 

KC=1. 4; 

RAg=-l-g/(s*ss); 

DA=DAO*(TP/TDAO)**RAG; 



% gas constant N-m/kg-K 
% Cp/Cv for air 
% calculation factor 
% AIR DENSITY AT PROJECTILE HEIGHT 



% Speed of Sound 
KR=KC*’'ss; 

% Vah=sqrt(KR-TP); % speed of sound at PI 



% curve fitting for Drag coefficient( experimental data input) 
x=( 0,0. 05,0. 15,0. 2,0. 23,0. 25,0. 3,0. 35,0. 4,0. 45,0. 5,0. 55,0. 6,. . . 

0. 65,0. 7,0. 75,0. 825,0. 875,0. 98,1. 0,1. 05,1. 1, 1. 2,1. 3,1. 4,1. 5,. . . 
1.75,2,2.5,2.9); 

\’=(0. 16,0. 16,0. 16,0. 16,0. 16,0. 16,0. 16,0. 16,0. 16 , 0. 16 , 0. 16 , 0. 16 , . . . 
b. 16,0. 16,0. 16,0. 16,0. 162,0. 17,0. 205,0. 225,0. 27,0. 282,0. 29,. . . 

0. 2875,0. 285,0. 275,0. 25,0. 23,0. 195,0. 17); 



cds=spline( x , y ) ; 



% spline curve fitting function 
% get spline coefficient 



% Trajectory Computation Begins 



t(l)=0; 

n=0; 

/o 

while Pln>a, n=n+l; 



% initial time 
main loop begin ***********: 






Update 

Pl=Pln; 

P2=P2n; 
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P3=P3n; 

P4=P4n; 

PVl=PVln; 

PAV2=PAV2n; 

PAV3=PAV3n; 

PV2=PAV2*P1; 

PV3=PAV3*P4; 

PVe=(PAV3-EAV)*P4; 

Height(n)=Pl -a; 

% velocity transformation matrix 

PR3=P3-EAV*t(n) ; % location relative to radar 

PG3=PR3-G3; % location relative to gun 

P2s=sin(P2); 

P2c=cos( P2); 

PR3s=sin(PR3); 

PR3c=cos( PR3) ; 

PG3s=sin(PG3); 

PG3c=cos( PG3) ; 

VTRll=P2s*PR3c; 

VTR21=P2s*PR3s; 

VTR31=P2c; 

VTR=(VTRll,P2c*PR3c, -PR3s; VTR21 ,P2c*PR3s ,PR3c; VTR31 , -P2s , 0) ; 

VTGll=P2s*PG3c; 

VTG21=P2s*PG3s; 

VTG31=P2c; 



% Relative Location 



PRXYZ=(P1*-VTR11; VTR21; VTR31)-(R4; 0; RZ); 
SLR=N0RM(PRXYZ); % SLANT RANGE FROM RADAR 

SRR(n)=SLR; 






■ 5 V ‘.V Vr ‘V -)V •k ■)> " 5 V ’>V -5V ^V ■5V Vc * “sV Vr -sV ■)> -.V Vf ‘V ;V Vr •> V ■>V •>- -.V Vr * ‘V Vc ' V ■5V Vr -sV •5V "sV * * •5V * V * 5 V -k k k k k k k k k 



DR(: ,n)=PRxvz/SLR; 

Dz=DR( 1 ,n)’^R2s+DR( 3 ,n)*R2c; 
Dn=DR(3,n)*R2s-DR(l,n)*R2c; 
De=DR( 2 ,n); 

traz(n)=atan2(Dn ,De)’’" 180/pi; 
trel(n ) =a s in (Dz)* 180/pi; 

' -V .V .V -V -V .V - V -V .t- -V .V -V 



% direction relative to radar 
% tracking direction, elevation 
% tracking direction, north comp. 

% tracking direction, east component 
% azimuthal angle 
?o elevation angle 

‘*-5V-5rVfVf'5VVrVfVcV?**VrVr‘;cV?***V*’5V****-5VVr’5V**‘.V*-5V-5V 



PGXYZ=DTG*(P1*-VTG11; VTG21; VTG31)-(G4; 0; GZ)); 
SL=norm(PGxyz); % slant range 

SRange(n)=SL; 

GRange(n)=a*acos((Pl**2+Gl**2-SL 2)/(2^’-Pl*Gl)); 

% ground range 



PGx=PGxyz( 1); 

PGz=PGxyz( 3); 

FRange( n)=a*atan( PGx/(Gl+PGz) ) ; % range in firing plane 

CRange(n)=PGxyz(2); % cross range from firing plane 



^Relative Velocity 
PV=(PV1; PV2; PV3); 

PVR=VTR‘>PV-(0; RV; 0); 

VR(: ,n)=PVR; % velocity relative to radar 

VRr(n)=PVl; 



% Time Interval and Doppler effect Estimation 
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VDopp=PVR'*PRxyz/SLR; 

Vdopl( n)=VDopp; 
fdl(n)=2*VDopp/ lambda; 
dtn=MB/VDopp; 
t(n+l)=t(n)+dtn; 

% Frictional Acceleration 

TP=TrefK+s*(Pl-Thl); 

Da=Dao*(TP/TDao)*-'-RAg; 

Vah=sqrt(KR*TP); 

PW=(PV1; PV2-WS; PVE-WE); 
PVUM=norm(PW); 

PWMM=PWM/Vah; 

PH=PW/PWM; 

HeadingR(: ,n)=VTR*PH; 

CD=ppva 1 ( cds , PWMM ) ; 

FA= -FAC*D a*CD*P WM*P W ; 

/o Location and Velocity Prediction 



% Doppler velocity 

% Doppler Frequency 
% square-law doppler detector 



% air temperature at the projectile 
% air density at projectile height 
% speed of sound at PI 
% PROJECTILE AIR VELOCITY 



% projectile heading 
?o heading as seen from radar 

% drag acceleration 



P A 1 = ( P V2** 2+P V3** 2 ) / P 1 + F A( 1 ) - g* ( a/ P 1 ) 2 ; 

PAA2=P2s*P2c*PAV3*-"2+FA(2)/Pl; 

PAA3=FA(3)/P4; 

PVln=PAl*dtn+PVl; 

Pln=0. 5-'(PVl+PVln)*dtn+Pl; 

PAV2n= ( PAA2--Mtn+PAV2 ) *( P 1 / P In )**2 ; 

P2n=0. 5*(PAV2+PAV2n)*dtn+P2; 

P4n=Pln‘"S in( P2n) ; 

PAV3n=( PAA3*dtn+PAV3)*( P4/P4n)**2; 

P3n=0. 5*(PAV3+PAV3n)*dtn+P3; 
end 



/o'* 

o 

/b 



*****•’;****** end of main loop ***-,v***5V****-v** 

Estimate Point of Impact 

Ih=Pl-a; 

IA=-PA1; 



IV=-PV1; 

dt=( sqrt( IV**2+2*Ih*IA) -IV)/IA; 
N=n+1; 

t(iN)=t(n)+dt; 

Pln=a; 

PVln=PAl*dt+PVl; 

PAV2n=(PAA2*dt+PAV2)*(Pl/Pln)**2; 
P2n=0. 5*(PAV2+PAV2n)*dt+P2; 
P4n=Pln*sin( P2n) ; 

PAV3n=( PAA3*dt+PAV3 )*( P4/P4n)**2; 
P3N=0. 5*( PAV3+PAV3N)*DT+P3; 



% Update 
Pl=Pln; 
P2=P2n; 
P3=P3n; 
P4=P4n: 
PVl=PVln; 
PAV2=PAV2n; 
PAV3=PAV3n; 
PV2=PAV2*P1; 
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PV3=PAV3*P4; 

Height(N)=0; 

% velocity transformation matrix 

PR3=P3-EAV*t(N); % location relative to radar 

PG3=PR3-G3; % location relative to gun 

P2s=sin(P2); 

P2c=cos( P2) ; 

PR3s=sin(PR3); 

PR3c=cos( PR3); 

PG3s=sin(PG3); 

PG3c=cos( PG3) ; 

VTRll=P2s*PR3c; 

VTR21=P2s*PR3s; 

VTR31=P2c; 

VTR=( VTRl 1 , P2C*PR3C , -PR3S; VTR2 1 , P2C*PR3S , PR3C; VTR31 , -P2S , 0) ; 

VTGll=P2s’'^PG3c; 

VTG21=P2s*PG3s; 

VTG31=P2c; 



% Relative Location 

PRXYZ=(Pl^-'-\TRll; VTR21; VTR31)-(-R4; 0; RZ)-; 
SRR(N)=norm(PRxyz); 



DR(: ,N)=PRxvz/SRR(N); 

Dz=DR( 1 ,N)"R2s+DR(3 ,N)*R2c; 
DN=DR ( 3 , N ) ’"^R 2 S -DR ( 1 , N ) *R 2 C ; 
De=DR(2,N); 

traz(N)=atan2(Dn ,De)“ 180/pi; 
trel( N)=asin( Dz)*l 80/pi; 



Radar parameters registration 



location relative to radar 
tracking direction, elevation 
TRACKING DIRECTION, NORTH COMP, 
tracking direction, east component 
azimuthal angle 
elevation angle 






PGxyz=DTG*(Pl*[ VTGll; VTG21; VTG31]-[GA; 0; Gz] ); 

SL=norm(PGxyz) ; % slant range 

SRange(N)=SL; 

GRange(N)=a*acos((Pl**2+Gl**2-SL 2)/( 2*P1*G1 ) ) ; 

7o ground range 

PGx=PGxvz( 1 ) ; 

PGz=PGxyz(3); 

FRange(N)=a*atan(PGx/(Gl+PGz)); % ground range in firing plane 
CRange(N)=PGxyz(2); % cross range from firing plane 



% Relative Velocity 
PV=[ PVl; PV2; PV3] ; 

PVR=VTR*PV-( 0; RV; 0]; 

VDopp=PVR' "PRxyz/SLR; % Doppler velocity 

Vdop 1 ( N ) =VDopp ; 

fdl(N)=2*VDopp/lambda; % Doppler Frequency 

VR(: ,N)=PVR; 

VRr(N)=PVl; 

PW=(PV1; PV2-WS; PVE-VE); 

HeadingRC : , N)=VTR*( PVT«7norm( PW,^ ) ; 
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for i=l:N; 

VRM( i)=norm( VR( : ,i)); 
VRrf(i)=VRr(i)/VRM(i); 

VRt( i)=( 1 -VRrf ( i ) )*VRM( i ) ; 
beamD( i) = l*pi/180’’^SRR( i); 
end 

EXETIMEl=etime( clock, t2) 

QUt pUt data 



•diary ww 
EXETIMEl 



MN=MN 



ng=N 

MB=MB 

tg=t(N) 

r=Pl 



dlatitude=(P2-G2)*180/pi 
dlongitude=PG3” 180/pi 
Rg=GRange(N) 
FRg=FRange(N) 
CRg=CRange(N) 

CRgmax=max( CRange) 
maxHeight=max( Height) 
maxfd=max( fdl) 
maxr ange=SRR( N) 
VDOPmin=min( Vdopl ) 
VDOPmax=max( Vdopl ) 
diary 



/o'' *' '' 






SAVE AND LOAD DATA FOR NEXT 



CALCULATION***** 



.'DEL 50010. MAT 

SAVE 50010 MAXRANGE EXETIMEl DR HEADINGR SRR HEIGHT T N FDl MB 
! DEL X50010, MAT 

SAVE X50010 TRAZ TREE VDOPl VRM VRRF VRT VR VRR T SRR BEAMD 



/o '' '* '' '* '* '' '* '* '' '' '' '* '' '* '* '' '' '' *' *' 

%PR 

%PRTRAK 

%XX0 

%RCSA 

%RANGE3 

%TDATA 



Execute 

0 / 

/o 

% 

0 / 

/o 

O/ 

/o 

0 / 

/o 

% 



other program 









PRINT OUT TRAJECTORY 

TRACKING PARAMETERS & OUTPUT 

RCS & RANGE ANALYSIS PROGRAM 

RCS SUBFUNCTION 

POWER REQUIREME.NT 

PROPAGATION FACTOR CURVE FITTING 



! del GH. met 



trajecto x*y 



plot(GRange, Height) ,grid; 
title( ' PROJECTILE TRAJECTORY'); 
xlabelC ' ground range (m)'); 
ylabel( ' projectile height (m)'); 
meta GH 



%prtsc 

%pause 

!del GC. met 



plot(GRange, CRange) ,grid; 
t it le( 'PROJECTILE TRAJECTORY'); 
xlabelC ' ground range (m)'); 
ylabelC ' Cross range (m)'); 
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meta GC 
% pause 
% prtsc 

! del RH. met 

plot(SRR,Height , ) ,grid; 
title(* PROJECTILE TRAJECTORY^); 
xlabel( ^ Slant Range (m)^); 
ylabel( ^ projectile height (m)*); 
meta RH 
% prtsc 
pause 

% !del TH. met 
plot(t , Height) , grid; 
xlabel( ^ flight time (sec)*); 
ylabel( * projectile height (m)*); 
meta TH 
% pause 
%prtsc 

!del TSGFC.met 
plot(t jSRange, 
xlabel( * flight 
ylabel( * slant , 
meta TSGFC 
%pause 
prtsc 

!del SGFCT. met 
plot(SRange,t , 
xlabelC * s lant , 
ylabel( * flight 
meta SGFCT 
%pause 
%prtsc 

!del TVR. met 
plot(t ,VRr , *+ ’ ) ,grid; 
xlabel( * flight time (sec)*); 
ylabel( * Vertical Speed* ); 
meta TVR 
%pause 
%prtsc 
!del TVH. met 
plot(t,VRt , *+' ) ,grid; 
xlabel( * Flight time (sec)*); 
ylabel( *Horizontal Speed* ); 
meta TVH 
%pause 
prtsc 

! del FH. met 

plot(FRange, Height ) ,grid; 

xlabel( * ground range in firing plane (m)*); 
ylabel( * projectile height (m)^); 
meta FH 
%pause 
%prtsc 
! del FC. met 

plot(FRange,CRange) ,grid; 
xlabel( * Firing plane Range (m)*); 



t jGRange, tjFRange, t ,CRange) ,grid; 
time ( sec) * ); 

ground, firing plane & cross ranges 



GRange,t , FRange,t , CRange,t) ,grid; 
ground, firing plane cross ranges 
time (sec) * ); 



(m)^); 



(m) ' 
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ylabel( ’ Cross range (tn)*); 
meta FC 
%pause 
prtsc 

! del GV. met 
plot(GRange,VRM) ,grid; 
xlabel( * Ground Range (m)^); 
ylabel( * Projectile Speed (m/sec)* ); 
meta GV 
%pause 
%prtsc 
! del GRV. met 
plot(GRange, VRrf ) ,grid; 
xlabel( ^ Ground range (ra)'); 

ylabel( ^ Ratio of radial velocity to speed'); 
meta GRV 



%pause 

%prtsc 

® ^ V: */r * * "i f : V -,V -sV * ‘V * * -k k k k 



TRAKING PARAMETERS CALCULATION AND OUTPUT ****** 



%**************** Tracking angular velocity ************ 
avell=diff(trel). /diff(t); 
avazl=dif f( traz). /diff(t); 

^^/^kkkkkkkkkkkkkkkkkkk PLOT kkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



for n=l: N-1 
ta(n)=t(n); 
end 

!del AVAEt,met 
plot(ta,avell, ’+' ), grid; 

title( 'RADAR TRACKING ELEVATION ANGULAR SPEED ’ ); 
xlabelC ' In-flight time (sec)’); 

ylabel( ' Elevation Angular velocity (degree/sec)’); 
meta AVAEt 
% pause 
prtsc 

! del AVAZt. met 

plot ( ta , avazl ,'+’), grid; 

titleC 'RADAR TRACKING AZIMUTHAL ANGULAR SPEED'); 
xIabeIC ' In-flight time (sec)'); 
ylabel( ' Azimuthal angular speed (degree)'); 
meta AVAZt 
%pause 
prtsc 

^^^kkkkkkkkkkkkkkkkkkkkkkkk Radar Tracking Parameter output ************ 
! del tret. met 
plot(t ,trel) ,grid; 

title( 'RADAR TRACKING ELEVATION ANGLE degree'); 
xlabel( ' In-flight time (sec)'); 
ylabel( ' Elevation Angle (degree)'); 

•meta thel 
% pause 
% prtsc 

! del trazl. met 
plot(tjtraz) ,grid; 

title( 'RADAR TRACKING AZIMUTHAL PARAMETER'); 
xlabel( ' In-flight time (sec)'); 
ylabel( 'Azimuthal angle (degree)'); 
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meta trazl 
%pause 
% prtsc 

! del AVAE. met 
plot(avell, ) jgrid; 

title( 'RADAR TRACKING ELEVATION ANGULAR VELOCITY '); 
xlabeK 'MEASUREMENT BASE 

ylabelC ' Elevation Angular velocity (degree/sec)'); 
meta AVAE 
% pause 
% prtsc 

!del AVAZ.met 
plot(avazl , ' +' ) jgrid; 

title( 'RADAR TRACKING AZIMUTHAL ANGULAR VELOCITY'); 
XlabeK 'MEASUREMENT BASE'); 
ylabeK 'AZIMUTH ANGULAR (degree)'); 



meta AVAZ 
%pause 
%prtsc 

EXETIMEl=etime( clock ,t2 ) 

/o ' \ 4 \ ^ 4 % 4 \ 4 \ 4 4 % 4 \ 4 % 4 % 4 



graph ‘iV‘V‘5V’3V*'>V'3V’3V*‘jVVr*’sWf “3V*Vr*';V';Wf ■>V*^V**';V*‘5V‘5Wr*sV‘5V^V 



!del TVD.met 
plot(t . fdl) ,grid; 
title( ^DOPPLER FREQUENCY'); 
xlabel( ' Flight time (sec)'); 
ylabelC ' Doppler Frequency (Hz)') 
meta TVD 
% pause 
prtsc 

! del treVD. met 
plot(trel,Vdopl) ,grid; 
t it le( 'DOPPLER Velocity'); 

xlabel( ' Radar tracking elevation angle (degree)'); 
ylabelC ' Doppler Velocity(m) ’ ) 
meta trehVD 
%pause 
% prtsc 

! del trazVD, met 

plot( traZjVdopl) ,grid; 

title( 'DOPPLER Velocity'); 

xlabel( ' Radar Tracking azimuthal Angle (degree)'); 
ylabelC ' Doppler Velocity(m) ' ) 
meta trazVD 
%pause 
%prtsc 

! del RVDl. met 
plot( SRR , Vdopl) ,grid; 
title( 'DOPPLER Velocity'); 
xlabel( ' Slant Range(m)'); 
ylabelC ' Doppler Velocity(m) ' ) 
meta RVDl 
%pause 
prtsc 

! del tVDl, met 
plot( t , Vdopl) jgrid; 
title( 'DOPPLER Velocity'); 
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xlabel( ’ In-flight t ime( sec) ’ ) 
ylabel( ^ Doppler Velocity( m) ’ ) 
meta tVDl 
prtsc 

Radar Tracking Parameter output ************ 

! del tret. met 
plot( t , trel) ,grid; 

title( ’RADAR TRACKING ELEVATION ANGLE degree’); 
xlabel( ’ In-flight time (sec)’); 
ylabel( ’ Elevation Angle (degree)’); 
meta thel 
% pause 
prtsc 

! del trazl. met 
plot( t ,traz) jgrid; 

title( ’RADAR TRACKING AZIMUTHAL PARAMETER’); 
xlabel( ’ In-flight time (sec)’); 
ylabel( ’Azimuthal angle (degree) ' ); 
meta trazl 
%pause 
prtsc 

* •/.' Vr “/r * * * * *.V * * * * * V: * * * ** * * ^ V * * Vr * : V * * * * * * * * * * * * * * Vr * * * * * * “* * * * * * * * * * * * * * 

%******* pQs and Power requirement Estimation******* 
load 50010 



t3=clock 

Bvl=l. 2*max(fdl); 

Bv2=l. 2*22031; 

Bv3=l. 2*12721; 

o^******vr*^v*vr***** radar parameters 
Los s= 10**0. 24; 

KT=4. 021E-21; 

SN=10; 

Lamda=3E5/ 10525E6; 

Fn= 10**0. 2; 

G= 10**4. 43; 



% video Bandwidth for 1st filter 
% video Bandwidth for 2nd filter 
% video Bandwidth for 3rd filter 
calculation v? ******************** 
% total loss 2. 4dB 



**** 



?o KT=K*T 

% minimum signal to noise ratio lOdB 
% wavelength meter 10525MHz 
% noise figure 2dB 
% TX or RX antenna gain 44. 3dB 



cons t=( ( 4*pi ) **3*KT^^Fn*Los s*SN ) / ( ( G*Lamda )2 ) ; 

Vr * ■/,* * * * * * V.' *.’.* * * * *.’.* Vr * * * * ** V * ‘V * * * * * Vr * * Vr * *.'.* * * *.V * * * * * *. 



% Radar cross section &. power requirement estimations 
for n=l: N 

if max(Height)— Height(n) 
ttop=t(n); 

Ltop=SRR(n); 
end 
end 

%******************* switch 1 power requirement 



% Summit point data 
% time at summit 
% range at summit 



Bcon=Bvl*const; 
for n=l: 1092 
SRRa(n)=SRR(n); 
ta(n)=t(n); 
uDl=DR(: ,n); 
h=HeadingR( : ,n); 
ct=h’*uDl; 
rcs=rcsa( ct) ; 
ct l(n)=ct; 
rcsl(n)=rcs; 



% res function 
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theta 1 ( n)=acos ( ct)*180/pi; 

Ll=SRR(n); 

PtO=Ll**4*Bcon/rcs; 

Ptl(n)=PtO; % power requied in free space 

if Ll=Ltop % required power at summit 

ptop=PtO; 
end 

SWitch 2 pOW€T rSC^uiretll^nt 

Bcon=Bv2*const; 
for n=1093; 1630 
SRRa(n)=SRR(n); 
ta(n)=t(n); 
uDl=DR( : ,n); 
h=HeadingR(: ,n); 
ct=h’ *uDl; 



rcs=rcsa(ct); % res function 

ct l(n)=ct; 
rcsl(n)=rcs; 

thetal(n)=acos( ct)* 180/pi; 

Ll=SRR(n); 

P 1 0 = L 1 * * 4 * B c o n / r c s ; 

Ptl(n)=PtO; 
if Ll=Ltop 
ptop=PtO; 
end 

switchl=max( Pt 1 ) ; 
end % end for loop 

o/ 3^d filter 



7o power requied in free space 
% required power at summit 



Bcon=Bv3*const; 
for n=l631: 1898 
SRRa(n)=SRR(n); 
ta(n)=t(n); 
uDl=DR( : ,n); 
h=HeadingR( : ,n) ; 
ct=h’ *uDl; 



rcs=rcsa( ct); 
ct l(n)=ct; 
rcsl(n)=rcs; 

thetal(n)=acos( ct)*180/pi; 
Ll=SRR(n); 
PtO=Ll**4*Bcon/rcs; 
Ptl(n)=PtO; 
if Ll=Ltop 
ptop=PtO; 
end 

switch2=max( Pt 1 ) ; 
end % end for loop 



O/ 



% 






'***** output 



% res function 



% power requied in free space 
% required power at summit 









! del wvl 
diary wwl 



N=N 



f lighttime=max( t) 
stoptime=max( ta) 
stoprange=max( SRRa) 
stopfd=fdl( 1898) 
t imetop=ttop 



65 



rangetop=Ltop 
powertop=ptop 
EXETI ME 2=etime( clock, t3) 

T0TALTIME=EXETIME1 + EXETIME2 
diary 

%**vc*^vvr****:vvr******yr*** graph radar cross section -j^**************** 
! del rcsl 1. met 
plot( SRRa, rcsl , ^ + * ) , grid; 

ylabel( * Radar cross section square meters*); 
xlabelC * Slant range* ) ; 
meta rcsll 
?oprtsc 
% pause 
! del rcsth. met 
p lot ( theta 1 , res 1 , * + * ) , grid; 

y label( I Radar cross section sc^uare meters*); 
xlabelC ' Aspect angle (degree)*^); 
meta rcsth 
% pause 
%prtsc 

! del rasp. met 

plotC ta, thetal , ’ + * ) , grid; 

ylabelC ’ Radar cross section sc^uare meters’); 

xlabelC ’ Aspect angle Cdegree) ); 

meta rasp 

%pause 

?oprtsc 

! del rasl. met 

plotCctl , ) ,grid; 

xlabelC * in-flight time Csec)’); 

ylabelC * Aspect angle Cdegree)*); 

meta rasl 

%prtsc 

/opause 

*•>*.> ‘V^V s a V e V a T i ab 1 e ^^'’VVr'jV^V^V'sWr'sV^'jVi'r-V'VVc’Vf 'sV’jV'jV'sV'jV^V'jV^'jV'jV^VVr^Wr'V^^Wf 



%\ del PwSOOlO. mat 

%save PwSOOlO ptlmax ctl t fdl rcsl LI 

VrVf printout pOW e r 

! del powerl.met 
plotCSRRa,Ptl ,’+*), grid; 
titleC*RADAR POWER REQUIREMENT *); 
xlabelC ’ Slant range Cm) ' 
ylabelC * Power requirement CWatts) * ); 
meta powerl 



N thetal Ptl SRR N 
requirement 



prtsc 

%pause 

***-.v,v-.v**** radar cross section CALCULATION Vr************•,v**^v** 






clear 

LOAD 50010 
for n=l:N 
uDl=DRC : ,n) ; 
h=HeadingRC: ,n); 
ct=h ' ‘'^uDl; 

/^function [ res , ct] =RCSC ct ) 
% constant declaration 
b=0. 155/2; 



% radar to target unit vector 
%unit heading vector of projectile 
% antenna beam angle 



% shell base radius 
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1=0. 5236; 
i=0. 4144; 
ta=b/l; 
ta2=ta*ta; 
ca2=l/(l+ta2); 
ca=sqrt( ca2); 
cx=abs(ct); 
cx2=cx-cx; 



% shell cone(nuzzle) height 
% shell cylinder height 
% =atan(b/l) i. e. half shell cone angle 



sx2=( l-cx2); 
sx=sqrt( sx2); 
rcsc=2'*’' i’'b-sx; 
b2=b*b; 

rcsb=pi‘'^b2*cx; 
if cx>=ca 
rcs=rcsc+rcsb; 
else 

rt=ta*cx/sx; 

rt2=rt*rt; 

rcs2=b2*’'( asin( rt )+rt*sqrt( l-rt2) ) - cx; 
rcsl=((b-i)/(sx2))-(l-cx2/ca2)‘-*'^(3/2); 
rcs=rcs l+rcs2+rcsc+rcsb/2; 

end % radar cross section case 2 

rcs0(n)=rcs; 
ctl(n)=ct; 

theta 1( n)=acos ( ct)* 180/pi; 
end 

! del res 11 

% plot(SRR,rcsO, ’+' ) jgrid; 

%ylabel( ’ Radar cross section square meters*); 

?oXlabel( * Slant range* ); 

%meta res 11 
%prtsc 
%pause 
! del rcsth 

p lot ( theta 1 , rcsO , * + * ) , grid; 

ylabelC * Radar cross section sc^uare meters*); 
xlabelC ' Aspect angle (degree)*); 
meta rcsth 
%pause 
! del rasp 

plot ( t , thetal , * + ' ) ,grid; 

ylabel( * Radar cross section s(^uare meters*); 
xlabel( * Aspect angle (degree) ); 
meta rasp 
/©pause 
! del ras 1 

plot(ctl, *+* ) ,grid; 
xlabel( ' in- flight time (sec)*); 
ylabel( ' Aspect angle (degree)’); 
meta rasl 



0 nd of program *’>V'5V’jV')V’>Vyr*'5V*Vr’»V**5VVryr'3iV*'5V')V*^Vyc’’>V'5Vyr^V'3V*>V’jV'5V'>V'5V 
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APPENDIX B. CROSS-SECTION ANALYSIS 



The radar cross-section has been discussed in Chapter 4. The derivation of the RCS 
equations is given below. The radar main beam direction is pointing in the direction 
D. The projectile axis (or heading) is pointing along h. Accordingly, the coordinate 
system {x,y.z) can be chosen so that x = D and the y-z plane is the plane on which the 
physical cross-sectional area of the projectile is to be projected. Another coordinate 
system (.v', r'. z') can be chosen so that = h and the projectile is a body of revolution 

around the .v' axis. These two coordinate systems can be set up jointly as follows ; 



II 


(5.1) 


X 

11 


(5.2) 


1 Z) X /! 1 




A A 

y = z X D 


(5.3) 


A ^ 

.v' = h 


(5.4) 




(5.5) 


II 

bj 

X 


(5.6) 


" r 

= z X h 


(5.7) 


A AAA 

[(/i X D]h -D] 


(5.8) 


A 


1 Z) X 1 




II 

1 

1 

‘ 


(5.9) 


1 >'p - > R 1 





Note that the coordinate (x.y.z) are used in this Appendix only and are not related 
to other parts of the thesis. 

The suface of the projectile is given by : 



6S 



/' + 2'' =/(x') 



(^. 10 ) 



The RCS is the area enclosed by projecting of this surface on the y-z plane. 

A. CASE 1 

For \ D • h\> cos a. Figure 25 on page 38. the conoid is hidded in the shadow re- 
gion. The projected areas can be easily seen to be : 



2 

= 3.14159-^ (5.11) 

(7, = I 5./H (5.12) 

A 2 = lab (5.13) 

= T, I 5./H (5.14) 

A 2 = A, (5.15) 

(73 = .43 I D./j I (5.16) 



where a is the diameter of the projectile and b is the cylinder height. The total projection 
area in ; 

Oj = <7] -f (72 -f (73 (5.17) 

B. CASE 2 

When \ D^ h \ <cosa, the geometr}’ of the target is shown in Figure 26 on page 
39. The RCS value consists of four different parts: A^ . A^, A^. . Each of them is given 

below : 

For area A,: 

= cV^ (5. IS) 

= (5.19) 



where 

c = tan a (5.20) 
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(c^ - tan^0^j^).v'^ + 2xp' xan^Oa.p = 0 


(B.2\) 


-'Cc - , , 

( - c‘) 


(5.22) 


Xq c" tan 

}c 2 -j 

( tan - c‘) 


(5.23) 


where .r/. _v/ is the maximum point of the ellipse. 




y y 

x' ~ tan 


(5.24) 


/(■'(' -yl = J ' tan + x' 


(5.25) 


X = V/ 


(5.26) 


(f‘.v' -f J-' tan d^^r) 

(tan'<?,,,-cV'' 


(5.27) 


At a point which is /■>, , x' = L, then: 




, 1 -' tan + Lc^ = 0 


(5.28) 


or 




r' 

tan 6^,^ 


(5.29) 


— Lc^ 

So the coordinate at p, is (-L , ^ , 0 ) on the z' axis, 

tan 

4 , 2 

■■ /• + - c-p 

tan 


(5.30) 


or 




cL(tanV"^"-c')^^'' 
tan O^^p 


(5.31) 
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Hence 



J, . zL(-4^)''= 



tan‘$ 



asp 



t1 , ^ In , 4J/2 ^ ^ 

= L c{ tan d 4- c ) ' 



tan^0 



asp 



Let .V be the normal unit vector of surface A{. 



A' = 



D cos'a - h cos 0 



asp 



2 4 40 

sin cos a + sin"a cos~6^^p 



then the projection area for is: 



tan a( ian'6,,p + tan a)''“( lanB,,, - tan a.y'\D • A1 



a, = 



tan'P 



asp 



gL , , _ x3/2 

• 2n ■> > 

Sin e^,p 



2 

cos Cf. 



(5.32) 



(5.33) 



(5.34) 



(5.35) 



(5.36) 



A 2 is part of a semi-circle (Figure 26 on page 39). A straight forward intergration gives: 



A, = 



a^{ sin '( 



tan y. 

I (>asp I 



) + 



tan a 
I tan e^^p 




n 



tan'o: 

tan'^^j^ 



} 



(5.37) 
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= a^{ sin '( 



tan S' 

I tan e^,p I 



) + 



tan a 
I tan e^,p I 




0 



tan'g 
tan B^^p 



]\b.h 



(5.3S) 



Ai is a rectangle: 
A-^ = lab 

so 

C 3 = A^ \ D 'K h \ 



(5.39) 



(5.40) 
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where a is the radius of the projectile base, and b is the height of the cylindrical part of 
projectile. 

At is a semicircle, hence : 

/J 4 = 3.14159-^ (5.41) 

= Ai\b>h\ (5.42) 

The total cross-section a, is the sum of upper 4 projection areas. 

(7, = ( 7 , 4- (7 2 -I- <73 -t- (74 (5.43) 
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APPENDIX C. PROPAGATION FACTORS 



Projectile 


Location 




Ground Range Height 


10m 


to 






impact point 




(m) 


(m) 


(dB) 


-482. 0000 


932. 4000 


0 


-382. 0000 


741. 2000 


0. 6000 


-333. 0000 


647. 9000 


-0. 7000 


-234. 0000 


455. 2000 


-0. 3000 


-182. 0000 


355. 7000 


-0. 9000 


-152. 0000 


294. 8000 


1. 9000 


-148. 0000 


289. 2000 


2. 3000 


-145. 0000 


283. 7000 


1. 5000 


-142. 0000 


278. 1000 


-0. 7000 


-140. 0000 


272. 6000 


-3. 8000 


-137. 0000 


267. 0000 


-3. 0000 


-132. 0000 


261. 5000 


0 


-131. 0000 


256. 0000 


2. 3000 


-129. 0000 


250. 4000 


3. 2000 


-125. 0000 


244. 8000 


2. 2000 


-123. 0000 


239. 3000 


-0. 8000 


-120. 0000 


233. 7000 


-5. 0000 


-117. 0000 


228. 2000 


-4. 2000 


-114. 0000 


222. 6000 


-0. 1000 


-111. 0000 


217. 1000 


3. 0000 


-108. 0000 


211. 5000 


3. 7000 


-105. 0000 


206. 0000 


2. 7000 


-103. 0000 


200. 4000 


-0. 3000 


-100. 0000 


194. 9000 


-6. 0000 


-97. 0000 


189. 3000 


-5. 3000 


-94. 0000 


183. 8000 


-14. 7000 


-91. 0000 


178. 2000 


3. 4000 


-88. 0000 


172. 6000 


2. 2000 


-85. 0000 


167. 1000 


3. 5000 


-83. 0000 


161. 5000 


0. 5000 


-80. 0000 


156. 0000 


-6. 4000 


-77. 0000 


150. 4000 


4. 3000 


-74. 0000 


144. 8000 


0. 7000 


-71. 0000 


139. 3000 


3. 7000 


-68. 0000 


133. 7000 


4. 6000 


-66. 0000 


128. 1000 


3. 8000 


-63. 0000 


122. 5000 


1. 1000 


-60. 0000 


117. 0000 


-6. 7000 


-57. 0000 


114. 4000 


-13. 4000 


-54. 0000 


105. 9000 


-0. 2000 


-52. 0000 


100. 3000 


4. 0000 


-49. 0000 


94. 7000 


5. 0000 


-46. 0000 


89. 2000 


4. 5000 


-43. 0000 


83. 6000 


2. 8000 


-40. 0000 


78. 0000 


1. 4000 



Antenna Height 

15m 20m 30m 50m 



(dB) 


(dB) 


(dB) 


(dB) 


0 


0. 6000 


-0. 3000 


-0. 3000 


6000 


0. 6000 


-0. 3000 


-0. 3000 


4000 


0. 4000 


-0. 3000 


-0. 3000 


1000 


1. 2000 


1. 2000 


-0. 2000 


0000 


1. 7000 


0. 7000 


1. 7000 


7000 


1. 2000 


-2. 6000 


1. 2000 


5000 


-3. 5000 


1. 5000 


-3. 4000 


0000 


0. 3000 


0 


0. 4000 


3000 


-1. 9000 


0. 3000 


-0. 7000 


4000 


-1. 0000 


-1. 9000 


-2. 1000 


0000 


-3. 0000 


-3. 2000 


-0. 9000 


4000 


2. 3000 


2. 6000 


2. 3000 


9000 


1. 5000 


-2. 4000 


-0. 3000 


7000 


-4. 1000 


1. 3000 


-2. 1000 


1000 


0. 2000 


0. 2000 


2. 2000 


6000 


3. 1000 


-1. 4000 


-0. 7000 


9000 


-4. 0000 


-5. 0000 


1. 6000 


9000 


-4. 0000 


-5. 0000 


1. 6000 


0000 


2. 9000 


3. 2000 


0. 7000 


5000 


2. 7000 


-0. 3000 


-0. 3000 


7000 


-5. 3000 


0. 8000 


-0. 3000 


7000 


-0. 9000 


1. 8000 


2. 7000 


7000 


3. 8000 


-2. 4000 


-2. 4000 


2000 


0. 6000 


3. 0000 


-0. 5000 


3000 


-5. 3000 


-5. 0000 


3. 0000 


7000 


-12. 5000 


3. 0000 


-11. 0000 


4000 


3. 2000 


1. 4000 


-1. 7000 


2000 


-3. 5000 


-0. 9000 


3. 1000 


5000 


-1. 5000 


3. 5000 


-0. 5000 


5000 


4. 2000 


-7. 5000 


-4. 8000 


4000 


2. 3000 


4. 1000 


2. 8000 


1000 


-10. 0000 


-3. 0000 


3. 0000 


9000 


2. 3000 


1. 9000 


-6. 9000 


9000 


4. 4000 


3. 0000 


3. 0000 


4000 


-1. 4000 


-4. 4000 


3. 2000 


4000 


-3. 3000 


4. 4000 


-6. 4000 


1000 


4. 1000 


-4. 7000 


1. 9000 


1000 


3. 6000 


2. 5000 


4. 0000 


4000 


0. 2000 


4. 6000 


4. 6000 


2000 


0. 2000 


4. 1000 


2. 7000 


8000 


4. 7000 


4. 7000 


3. 1000 


0000 


2. 4000 


-3. 9000 


-6. 9000 


2000 


-14. 4000 


2. 5000 


3. 6000 


5000 


1. 5000 


3. 5000 


1. 5000 


2000 


4. 7000 


-10. 2000 


-2. 2000 



0 

- 0 . 

- 1 , 

1 . 

-3, 

0 . 

1 

1 , 

-2 

-3, 

1 . 

2 

0 , 

-5. 

- 1 , 

2 . 

2 , 

- 1 , 

- 6 , 

0 . 

3 

2 , 

- 2 , 

-5. 

-4, 

4 

4, 

-5, 

-3, 

3 

7 

0 

-9 

-0 

4, 

4, 

0 

-4 

0 , 

4 

4 

0 

14 

0 
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- 37. 0000 
- 34. 0000 
- 31 . 0000 
- 29. 0000 
- 26. 0000 
- 23. 0000 
- 20. 0000 
- 17. 0000 
- 14. 0000 
- 11. 0000 
- 9. 0000 
- 6. 0000 
- 3 . 0000 
0 



72. 5000 

66. 9000 
61 . 3000 

55. 7000 
50. 2000 

44. 6000 
39. 0000 

33. 4000 

27. 9000 

22. 3000 

16. 7000 

11. 1000 

5. 5000 
0 



- 12 . 4000 
- 4 . 8000 

1 . 4000 

3. 7000 

4. 7000 

3. 5000 

2. 7000 
1 , 8000 

1. 1000 

- 0. 1000 
1 . 0000 
- 0. 6000 

4. 7000 
- 2. 6000 



4 . 2000 

4. 7000 

1 . 4000 
- 7. 0000 
- 4. 6000 

2. 3000 
3 . 7000 

2 . 4000 
2 . 4000 

0. 7000 

- 0. 3000 
- 0. 9000 

3. 6000 
- 4. 8000 



16. 0000 
- 12. 6000 

1 . 4000 

4. 3000 

2. 3000 
- 6.4000 
- 2. 1000 

2. 4000 

3 . 7000 
2. 2000 
0. 3000 
- 0. 2000 

2 . 7000 
- 5. 0000 



4. 6000 
0 . 9000 

- 3. 1000 

4. 5000 
- 1. 9000 
- 1. 4000 

4. 4000 
0 

- 1. 7000 

3. 6000 

3. 3000 

- 0. 2000 

2. 4000 
- 4. 6000 



4. 6000 
- 3 . 8000 

2. 9000 

2 . 6000 
- 3. 3000 

4. 5000 
- 11 . 1000 

4. 5000 
- 9. 0000 
- 10. 7000 

4. 1000 

4. 1000 
2. 8000 
- 3. 2000 
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APPENDIX D. MISSED DETECTION AND FALSE ALARM 



clear 

MN=1000; 

Npeak=MN/2; 
dtmax=0. 08650; 

Bvl=63900; 

FNnoise=Bvl*dtraax/Npeak; 

SNR0=10. 0; 
s tep=0. 5; 

Nf=5; 

Mf=29; 

thshstep=0, 005; 
thsh0=0. 355; 
for m=l:Mf 

thsh=m“ thshstep+thshO; 
thresh(m)=thsh; 

FNEnoise(m)=FNnoise*2*asin( sqrt ( thsh) ) /pi; 
end 

x=SNR0-step; 
for n=l: Nf 
x=x+step; 

SNR(n)=x; 
sgnr=10**(0. l*x); 

% transfer the signal-to-noise parameter to be used in fdetect.m 
snr=sgnr 

! del hm tempsnr. mat 
save hm tempsnr snr 

for m=l: Mf 
thsh=thresh(m); 
pfan=exp( -thsh*sgnr); 

Pfa(n,m)=pfan; 

FNerrorC n,m)=pf an*FNEnoise(m); 
end 

Pmiss(n , l)=Pmissdct ( sgnr , 0 , thresh( 1) ) ; 
for m=2:Mf 

Pmissd=Pmissdct( sgnr , thresh( m- 1) , thresh(m) ) ; 

Pmiss(n ,m)=Pmiss(n,m- 1)+Pmissd; 
end 

end 

! del hm pfaulse. mat 

save hm pfaulse Nf Mf thresh FNEnoise SNR Pmiss FNerror Pfa 

! del hm pfaulse 
diary hm pfaulse 

disp( * Threshold Missed Noise Faulse') 

disp ( ' Detection Error Alarm ^ ) 

for n=l:Nf 
SN_rat io=SNR(n) 

(thresh; Pmiss(n,:); FNerrorCn,: ); Pfa(n,:))^; 
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end 

diary 

pause 

!del hm pfaulse.met 

plot (thresh, Pmiss( 1 ,:),* + \ thresh, FNerror( 1 , : ) , * o * ) , grid; 

tit le( * Signal-to-Noise Ratio = 10.0 db* ) 

xlabel( * Threshold-to-Signal Power Ratio*); 

ylabel( * Missed Detection ( + ) &. Noise Error (o)*) 

meta hm pfaulse 

pause 

plot (thresh, Pmiss( 2 , thresh, FNerror( 2 , : ) , * o* ) , grid; 

tit le( * Signal-to-Noise Ratio = 10.5 db* ) 

xlabel( *Threshold-to-Signal Power Ratio'); 

ylabel( ' Missed Detection (+) & Noise Error (o)') 

meta hm pfaulse 

pause 

plot (thresh, Pmiss( 3 ,:),’+' , thresh, FNerror ( 3 , : ) , ' o ' ) , grid; 

title( ' Signal-to-Noise Ratio = 11.0 db' ) 

xlabel( ' Threshold- to-Signal Power Ratio*); 

ylabel( ’ Missed Detection ( + ) &. Noise Error (o)*) 

meta hm pfaulse 

pause 

plot (thresh, Pmiss( 4 ,:),*+* , thresh, FNerror( 4 , : ) , * o * ) , grid; 

tit le( * Signal-to-Noise Ratio = 11.5 db*) 

xlabel( * Threshold-to-Signal Power Ratio*); 

ylabel( * Missed Detection ( + ) &. Noise Error (o)*) 

meta hm pfaulse 

pause 

plot (thresh, Pmiss( 5 ,:),*+* , thresh, FNerror(5 , : ) , * o* ) ,grid; 

tit le( * Signal -to-Noise Ratio = 12.0 db* ) 

xlabel( * Threshold-to-Signal Power Ratio'); 

ylabel( ’ Missed Detection ( + ) &. Noise Error (o)*) 

meta hm pfaulse 

pause 






function pmsd = Pmissdct ( snr , thsha, thshb) 

% The parameter snr for the singal-to-noise ratio (NOT in db) should 
% have been saved to hm tempsnr. mat. 

%!del hm tempsnr mat 
%save hm tempsnr snr 

pmsd=quad( * f detect ’ , thsha , thshb, 0. 0001); 






Vr Vr Vr tV tW? V* Vr Vr Vr Vr Vr Vc Vc V? vWt Vr Vr Vc Vr **V.'*************Vr** 



% The parameter sgnr has to be specified in the program, 
function pdens = fdetect(v) 
load hm tempsnr 
ziv=2*snr*sqrt( v)*i; 

LIoz=log(BESSELN(0,ziv)); 

Lpdens=LIoz+log( snr ) -snr*( v+1) ; 
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pdens=exp( Lpdens ) ; 

end of program ^v*-)V****-**********^v*****-5V**'3V*‘5V-)V*jv*** 

^'jV'/V'jVVf'jV'sV'jV'sV’sVVc'jV'jVVf^Vc^V'sV OU t pU t data **'*^'f **VVc'»VVf’5VV“'5V*'5V’;V****';V’sV'5V'jV**VVf'5V'sV*'sV’>V'3V’>V*'3V*V';V* 



Threshold Missed Noise Faulse 

Detection Error Alarm 

S/N = 10. 0 db 



0. 3600 


0. 0266 


0. 1237 


0. 0273 


0. 3650 


0. 0278 


0. 1187 


0. 0260 


0. 3700 


0. 0290 


0. 1138 


0. 0247 


0. 3750 


0. 0303 


0. 1091 


0. 0235 


0. 3800 


0. 0316 


0. 1046 


0. 0224 


0. 3850 


0. 0329 


0. 1002 


0. 0213 


0. 3900 


0. 0343 


0. 0961 


0. 0202 


0. 3950 


0. 0357 


0. 0921 


0. 0193 


0. 4000 


0. 0371 


0. 0883 


0. 0183 


0. 4050 


0. 0386 


0. 0846 


0. 0174 


0. 4100 


0. 0401 


0. 0810 


0. 0166 


0. 4150 


0. 0417 


0. 0777 


0. 0158 


0. 4200 


0. 0433 


0. 0744 


0. 0150 


0. 4250 


0. 0449 


0. 0713 


0. 0143 


0. 4300 


0. 0466 


0. 0683 


0. 0136 


0. 4350 


0. 0483 


0. 0654 


0. 0129 


0. 4400 


0. 0501 


0. 0627 


0. 0123 


0. 4450 


0. 0519 


0. 0600 


0. 0117 


0. 4500 


0. 0537 


0. 0575 


0. 0111 


0. 4550 


0. 0556 


0. 0551 


0. 0106 


0. 4600 


0. 0575 


0. 0527 


0. 0101 


0. 4650 


0. 0595 


0. 0505 


0. 0096 


0. 4700 


0. 0615 


0. 0484 


0. 0091 


0. 4750 


0. 0636 


0. 0463 


0. 0087 


0. 4800 


0. 0656 


0. 0443 


0. 0082 


0. 4850 


0. 0678 


0. 0424 


0. 0078 


0. 4900 


0. 0699 


0. 0406 


0. 0074 


0. 4950 


0. 0721 


0. 0389 


0. 0071 


0. 5000 


0. 0744 


0. 0372 


0. 0067 



S/N = 10. 5 db 



0. 3600 


0. 0211 


0. 0798 


0. 0176 


0. 3650 


0. 0221 


0. 0760 


0. 0166 


0. 3700 


0. 0232 


0. 0724 


0. 0157 


0. 3750 


0. 0243 


0. 0690 


0. 0149 


0. 3800 


0. 0255 


0. 0658 


0. 0141 


0. 3850 


0. 0267 


0. 0627 


0. 0133 


0. 3900 


0. 0279 


0. 0597 


0. 0126 


0. 3950 


0. 0291 


0. 0569 


0. 0119 


0. 4000 


0. 0304 


0. 0542 


0. 0112 


0. 4050 


0. 0318 


0. 0516 


0. 0106 


0. 4100 


0. 0332 


0. 0491 


0. 0100 


0. 4150 


0. 0346 


0. 0468 


0. 0095 


0. 4200 


0. 0360 


0. 0446 


0. 0090 


0. 4250 


0. 0375 


0. 0424 


0. 0085 


0. 4300 


0. 0391 


0. 0404 


0. 0080 


0. 4350 


0. 0406 


0. 0385 


0. 0076 
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0 . 4400 




0. 0423 


0 . 0366 


0. 0072 


0 . 4450 




0. 0439 


0. 0349 


0. 0068 


0 . 4500 




0 . 0456 


0 . 0332 


0. 0064 


0. 4550 




0. 0474 


0. 0316 


0 . 0061 


0 . 4600 




0. 0492 


0. 0301 


0. 0057 


0 . 4650 




0. 0510 


0. 0286 


0 . 0054 


0 . 4700 




0. 0529 


0. 0272 


0 . 0051 


0. 4750 




0. 0548 


0. 0259 


0. 0048 


0 . 4800 




0. 0567 


0. 0247 


0. 0046 


0. 4850 




0. 0588 


0. 0235 


0 . 0043 


0. 4900 




0. 0608 


0. 0223 


0. 0041 


0. 4950 




0. 0629 


0. 0213 


0 . 0039 


0. 5000 




0. 0650 


0 . 0202 


0. 0037 


= 11 . 0 


db 








0 . 3600 




0 . 0163 


0. 0487 


0. 0108 


0 . 3650 




0. 0172 


0. 0461 


0. 0101 


0 . 3700 




0 . 0181 


0. 0436 


0. 0095 


0 . 3750 




0 . 0191 


0. 0413 


0. 0089 


0 . 3800 




0. 0201 


0 . 0391 


0. 0084 


0. 3850 




0. 0211 


0. 0370 


0. 0079 


0. 3900 




0. 0222 


0 . 0350 


0. 0074 


0. 3950 




0. 0233 


0 . 0331 


0. 0069 


0.4000 




0. 0244 


0. 0313 


0. 0065 


0 . 4050 




0. 0256 


0. 0296 


0 . 0061 


0. 4100 




0. 0268 


0. 0280 


0. 0057 


0. 4150 




0. 0281 


0. 0265 


0. 0054 


0. 4200 




0. 0294 


0 . 0251 


0 . 0051 


0. 4250 




0. 0307 


0. 0237 


0. 0047 


0. 4300 




0. 0321 


0. 0224 


0. 0045 


0 . 4350 




0. 0335 


0. 0212 


0. 0042 


0 . 4400 




0. 0350 


0 . 0201 


0. 0039 


0. 4450 




0. 0365 


0 . 0190 


0. 0037 


0 . 4500 




0. 0380 


0 . 0179 


0. 0035 


0. 4550 




0. 0396 


0 . 0170 


0. 0033 


0. 4600 




0. 0413 


0. 0160 


0 . 0031 


0. 4650 




0. 0430 


0 . 0151 


0. 0029 


0 . 4700 




0. 0447 


0 . 0143 


0. 0027 


0 . 4750 




0. 0465 


0. 0135 


0. 0025 


0. 4800 




0. 0483 


0 . 0128 


0. 0024 


0 . 4850 




0. 0502 


0 . 0121 


0. 0022 


0 . 4900 




0 . 0521 


0. 0114 


0 . 0021 


0 . 4950 




0. 0540 


0 . 0108 


0. 0020 


0 . 5000 




0. 0561 


0 . 0102 


0. 0018 


= 11.5 


db 








0. 3600 




0. 0123 


0 . 0280 


0 . 0062 


0. 3650 




0. 0130 


0. 0263 


0. 0058 


0 . 3700 




0. 0138 


0 . 0247 


0 . 0054 


0. 3750 




0. 0146 


0. 0232 


0. 0050 


0 . 3800 




0 . 0154 


0. 0218 


0 . 0047 


0. 3850 




0. 0163 


0. 0205 


0 . 0043 


0. 3900 




0 . 0172 


0 . 0192 


0 . 0041 


0. 3950 




0. 0182 


0 . 0181 


0 . 0038 
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0. 4000 


0 . 0191 


0. 0169 


0. 0035 


0. 4050 


0. 0201 


0. 0159 


0. 0033 


0 . 4100 


0 . 0212 


0 . 0149 


0. 0031 


0 . 4150 


0 . 0223 


0. 0140 


0. 0028 


0 . 4200 


0 . 0234 


0. 0132 


0. 0027 


0.4250 


0. 0246 


0. 0123 


0. 0025 


0. 4300 


0 . 0258 


0 . 0116 


0. 0023 


0. 4350 


0 . 0271 


0. 0109 


0 . 0021 


0. 4400 


0 . 0284 


0. 0102 


0. 0020 


0. 4450 


0. 0297 


0. 0096 


0 . 0019 


0 . 4500 


0 . 0311 


0. 0090 


0 . 0017 


0 . 4550 


0. 0325 


0. 0084 


0. 0016 


0. 4600 


0. 0340 


0. 0079 


0. 0015 


0 . 4650 


0 . 0355 


0 . 0074 


0. 0014 


0 . 4700 


0 . 0371 


0. 0070 


0. 0013 


0.4750 


0. 0387 


0. 0065 


0. 0012 


0 . 4800 


0. 0404 


0 . 0061 


0. 0011 


0. 4850 


0 . 0421 


0 . 0057 


0. 0011 


0. 4900 


0 . 0439 


0. 0054 


0. 0010 


0.4950 


0. 0457 


0. 0050 


0. 0009 


0 . 5000 


0. 0476 


0 . 0047 


0. 0009 


S/N = 12 . 0 db 


0. 3600 


0. 0090 


0 . 0151 


0. 0033 


0 . 3650 


0 . 0096 


0. 0140 


0. 0031 


0 . 3700 


0 . 0102 


0 . 0131 


0. 0028 


0 . 3750 


0 . 0108 


0 . 0122 


0. 0026 


0. 3800 


0. 0115 


0. 0113 


0. 0024 


0 . 3850 


0 . 0122 


0 . 0105 


0. 0022 


0. 3900 


0. 0130 


0. 0098 


0 . 0021 


0. 3950 


0 . 0138 


0 . 0091 


0. 0019 


0. 4000 


0 . 0146 


0. 0085 


0. 0018 


0. 4050 


0 . 0154 


0 . 0079 


0. 0016 


0 . 4100 


0 . 0163 


0 . 0074 


0. 0015 


0. 4150 


0. 0173 


0. 0069 


0. 0014 


0. 4200 


0 . 0182 


0. 0064 


0. 0013 


0. 4250 


0. 0192 


0. 0059 


0 . 0012 


0 . 4300 


0. 0203 


0. 0055 


0. 0011 


0. 4350 


0 . 0214 


0. 0051 


0. 0010 


0 . 4400 


0. 0225 


0 . 0048 


0. 0009 


0 . 4450 


0 . 0237 


0 . 0044 


0. 0009 


0. 4500 


0 . 0249 


0. 0041 


0. 0008 


0. 4550 


0 . 0261 


0. 0038 


0. 0007 


0. 4600 


0. 0275 


0. 0036 


0. 0007 


0. 4650 


0 . 0288 


0 . 0033 


0. 0006 


0 . 4700 


0 . 0302 


0. 0031 


0. 0006 


0. 4750 


0 . 0317 


0. 0029 


0. 0005 


0. 4800 


0. 0331 


0. 0027 


0. 0005 


0. 4850 


0. 0347 


0. 0025 


0. 0005 


0 . 4900 


0. 0363 


0. 0023 


0. 0004 


0. 4950 


0. 0379 


0. 0022 


0. 0004 


0. 5000 


0. 0396 


0. 0020 


0. 0004 


VfsV Vr "jV* Vt* "jV* ■jV'jV 


end of data 
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